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ABSTRACT 

We analyze data from the Hubble Space Telescope's Advanced Camera for Surveys of the globular 
cluster Omega Cen. We construct a photometric catalog of 1.2 x 10® stars over a 10' x 10' central field 
down to below i?F435W = 23 (M ~ 0.35Mq). The 2.5- to 4-year baseline between observations yields 
a catalog of some 10^ proper motions over a smaller area, with 53,382 "high-quality" measurements in 
a central i? < 2' field. Artificial-star tests characterize the photometric incompleteness. We determine 
the cluster center to ~ 1" accuracy from star counts using two different methods, one based on iso- 
density contours and the other on "pie-slices" . We independently confirm the result by determining 
also the kinematical center of the HST proper motions, as well as the center of unresolved light seen 
in 2MASS data. All results agree to within their 1-2" levels of uncertainty. The proper-motion 
dispersion of the cluster increases gradually inwards, but there is no variation in kinematics with 
position within the central ~ 15": there is no dispersion cusp and no stars with unusually high 
velocities. We measure for the first time in any globular cluster the variation in internal kinematics 
along the main sequence. The variation of proper-motion dispersion with mass shows that the cluster 
is not yet in equipartition. There are no differences in proper-motion kinematics between the different 
stellar populations of Omega Cen. Our results do not confirm the arguments put forward by Noyola, 
Gebhardt & Bergmann to suspect an intermediate-mass black hole (IMBH) in Omega Cen. They 
determined linc-of-sight velocity dispersions in two 5" x 5" fields, and reported higher motions in their 
central field. We find the proper-motion kinematics to be the same in both fields. Also, we find that 
they (as well as other previous studies) did not accurately identify the cluster center, so that both 
of their fields are in fact 12" from the true center. We also do not confirm the central density cusp 
they reported (in part due to the different center, and in part due to biases induced by their use of 
unresolved light). The surface number-density distribution near the center does not differ strongly 
from a single-mass King model, although a shallow cusp may not be ruled out. In Paper II we present 
new dynamical models for the high-quality data presented here, with the aim of putting quantitative 
constraints on the mass of any possible IMBH. 

Subject headings: catalogs — globular clusters: general — globular clusters: individual (Omega Cen) 
— stars: kinematics — techniques: image processing, photometric. 



1. INTRODUCTION 

There is much current interest in finding black holes 
at the centers of globular clusters (GCs). Supermassive 
(> IQ^Mq) black holes are known to exist at the cen- 
ters of galaxies, and it has been demonstrated that the 
black-hole mass is correlated with the velocity dispersion 
of the galaxy. If GCs house massive central black holes 
that follow the same correlation, then based on their cen- 
tral velocity dispersions of order '--^lO km/s, we would 
expect to find an intermediate-mass black hole (IMBH) 
with a mass of a few thousand solar masses (Gebhardt, 
Rich & Ho 2002). Plausible formation scenarios exist to 
form such IMBHs in the centers of globular clusters from 
realistic initial conditions (Portegies Zwart & McMillan 
2002). Determining whether globular clusters harbor 
IMBHs will answer important questions about how clus- 
ters form, and about what kinds of circumstances give 
rise to IMBHs (e.g., van der Marel 2004). 

'BASED ON OBSERVATIONS WITH THE NASA/ESA HUB- 
BLE SPACE TELESCOPE, OBTAINED AT THE SPACE TELE- 
SCOPE SCIENCE INSTITUTE, WHICH IS OPERATED BY 
AURA, INC., UNDER NASA CONTRACT NAS 5-26555. 



An IMBH should induce a power-law cusp in the 
stellar-density profile (Baumgardt, Makino & Hut 2005). 
Such cusps are not uncommon in globular clusters (Noy- 
ola & Gebhardt 2006). However, various stages of core 
collapse can introduce similar cusps, so observing a cusp 
is not definitive evidence for an IMBH. X-ray or radio 
emission may point towards accretion onto an IMBH 
(Pooley & Rappaport 2006; Kong 2007; Ulvestad, Greene 
& Ho 2007), but other explanations for the emission are 
difficult to rule out. Moreover, globular clusters generally 
have little gas to accrete. The most unambiguous way 
to identify and weigh an IMBH is therefore to find the 
signature it induces in the kinematics of nearby stars. 
There arc two ways in which kinematics could identify 
an IMBH. First, there may be a general increase in the 
velocity dispersion towards the center that cannot be ac- 
counted for by the visible matter. Second, there may be 
stars moving faster than would be allowed by the clus- 
ter's nominal escape velocity (Drukier & Bailyn 2003). 
One might even hope to observe stars in Keplerian orbit, 
as has been seen at the center of our own Galaxy (Genzel 
et al. 2003; Chez et al. 2005). 

Kinematical evidence for IMBHs on the basis of an in- 
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crease in the line-of-sight velocity dispersion towards the 
center has been presented for the globular clusters M15 
and Gl (van der Marel et al. 2002; Gerssen et al. 2002; 
Gebhardt, Rich, & Ho 2002, 2005). However, in M15 the 
implied dark mass could also be attributed to segregation 
of dark remnants towards the center (e.g., Baumgardt 
et al. 2003a), while in Gl the statistical significance of 
the implied dark mass is not strong (Baumgardt et al. 
2003b). The limitations of these past investigations are 
due in part to their use of line-of-sight velocities. Proper 
motions can generally put more powerful constraints on 
an IMBH than can line-of-sight velocities, because they 
can be measured for many more stars. Since proper mo- 
tions probe two components of the motion, they also 
yield better constraints on the velocity anisotropy and 
therefore help break the well-known degeneracy between 
mass and anisotropy. In cases where one can assume 
spherical symmetry, this degeneracy is resolved com- 
pletely (Leonard & Merritt 1989). Furthermore, while 
spectroscopic line-of-sight velocity studies are limited to 
the bright stars, the Hubble Space Telescope (HST) can 
observe proper motions for the much more plentiful main- 
sequence stars. This gives better statistics to probe closer 
to the center. IMBH limits from HST proper motions 
have been presented for M15 (McNamara et al. 2003 and 
van den Bosch et al. 2006) and 47 Tuc (Mcaughlin et 
al. 2006). These studies relied, at least in part, on data 
from the Wide Field and Planetary Camera 2 (WFPC2) 
instrument. Several studies have since been started or 
planned to study this problem exclusively with data from 
newer Advanced Camera for Surveys (ACS) or Wide 
Field Camera 3 (WFC3) instruments for larger samples 
of clusters (e.g., GO-9835 and GO-10474, PI-Drukier; 
GO-10401/10841, PTChandar; GTO-10335/11801, PI- 
Ford; and GO-11609, Pl-Chaname). Two epochs are al- 
ready in-hand for several of the target clusters in these 
samples and results may start to come out soon. 

Recently, Noyola, Gebhardt, and Bergmann (2008, 
NGB08) presented a new study of the globular clus- 
ter Omega Cen. They used Gemini IFU spectroscopy 
to measure the line-of-sight velocity dispersion of unre- 
solved light in two 5" x 5" fields, one at the cluster cen- 
ter and one at i? = 14" from the center. The dispersion 
in the central field (23.0 ± 2.0 km/s) exceeded that in 
the off-center field (18.6 ± 1.6 km/s). Based on this in- 
crease towards the center they argued for the presence of 
an IMBH of mass 4.0t?;^^ x IO^'Mq. NGB08 also mea- 
sured the surface-brightness profile of unresolved light 
from HST/ ACS images. They found it to have a shallow 
central cusp of logarithmic slope 7 = 0.08 ± 0.03, also 
consistent with the presence of an IMBH. 

A considerable amount of ACS data of Omega Cen al- 
ready exists in the HST Data Archive. In the present 
paper we collect and analyze these data to obtain pho- 
tometric and proper-motion catalogs for very large num- 
bers of stars (of order 10^ photometric measurements and 
10^ proper motions, respectively). We use these catalogs 
to perform detailed star-count and kinematical analyses. 
The present paper is confined to the observational do- 
main as much as possible, and we focus on testing and 
augmenting the results of NGB08. Detailed dynamical 
modeling of the new data is presented in a companion 
paper (van der Marel & Anderson 2009; Paper II). 

This paper is organized as follows. In Section[2]we give 



TABLE 1 

Data available for the proper-motion study. 



Dataset 


CENTRAL FIELD 


OUTER FIELD 


GO-9442 
(2002.49) 


3 X 340s F435W 

3 X 340s F625W 

4 X 440s F658N 
(Pointing 5) 


3 X 340s F435W 

3 X 340s F625W 

4 X 440s F658N 
(Pointings 1 and 4) 


GO-10252 
(2004.95) 




5 X 340s F606W 
5 X 340s F658N 


GO-10755 
(2006.56) 


4x80s F606W 
4x90s F814W 





an overview of the three ACS/WFC data sets that are 
part of our study. In Section [3] we construct a photo- 
metric catalog covering the inner part of the cluster, and 
we determine proper motions for the stars that could 
be measured well in two epochs. In Section |4] we deter- 
mine the position of the cluster center using our new star 
lists and proper motions, as well as ground-based 2MASS 
data. In Section [S] we derive the number-density profile. 
In Section[n]we do an initial analysis of the proper-motion 
kinematics (with more detailed analysis following in Pa- 
per II), and we calculate the proper- motion dispersions 
in the fields observed by NGB08. Our results for the 
cluster center, density profile, and kinematical gradient 
are all different from those presented by NGB08. We 
do not confirm the arguments put forward by them to 
suspect an IMBH in Omega Cen. In Section [7] we com- 
bine our proper motions and photometry to determine 
whether the cluster core is in equipartition and to de- 
termine whether the different stellar populations in the 
cluster have different kinematics. Finally, in Section [H] 
we present the conclusions of our study. 

2. OBSERVATIONS 

This paper draws on three HST data sets, listed in Ta- 
ble [H All three were taken with the Wide-Field Chan- 
nel (WFC) of the Advanced Camera for Surveys (ACS), 
which is made up of two 2048 x 4096-pixel detectors with 
a pixel scale of about 50 mas / pixel, covering a rhombus- 
shaped field that is roughly 3!4 on a side. 

The first data set is GO-9442 (Pl-Cool). It was taken 
over the period between 27 and 30 Jun 2002 and consists 
of a 3x3 mosaic of pointings, centered on the cluster 
center. Each pointing has 3 deep exposures and 1 shal- 
low exposure in each of F435W (i?F435w) and F625W 
(^F625w) and 4 deep observations through F658N (H- 
alpha, -ffpessN)- Each of the deep observations was offset 
relative to the others, so that no star would fall in the 50- 
pixel-tall inter-chip gap in more than one deep exposure 
for each filter. The mosaic of pointings covers the inner 
lO'xlO', extending out to about two core radii (fc ~ 2!5). 
The left panel of Figure [T] shows the field coverage for the 
-Bf435w filter. HST mis-pointed in the last pointing (up- 
per right), and there is a small gap in the field coverage 
there. We were careful to include the effect of this gap 
in our star-count analysis. 

The second data set is GO-10252 (PI- Anderson). This 
program was taken on 11 Dec 2004, and consists of 5 deep 
exposures and one shallow exposure in each of F606W 
(Vf606w) and F814W (/f814w)- The goal of these ob- 
servations was to construct accurate PSF models and an 
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Fig. 1. — At each point in the field centered on Omega Cen, 
we show the number of exposures available for each of the two 
broadband colors for the three data sets. The circles indicate one 
core radius (2^5) and two core radii. See text for a description 
of the data in each program. The baseline between the first and 
second panel is 2.46 years, and between the first and third is 4.07 
years. 

accurate distortion solution for the two most commonly 
used filters, thus the various exposures have large offsets 
of about 500 pixels between them. Unfortunately, at the 
time these observations were taken, the camera was con- 
siderably out of focus (it was adjusted soon afterwards), 
and the images were therefore not useful for their original 
aim of probing the "typical" PSF behavior and distortion 
solution. However, they are still useful for our purposes 
here, since we can tailor-make PSFs for these images, and 
can measure our motions locally, so that small errors in 
the global distortion solution will not impact our results. 
The middle panel of Figured] shows that this field is cen- 
tered near the south-east corner of the mosaic, along the 
cluster's major axis. The baseline between this and the 
first dataset is 2.49 years. 

The third data set is GO-10775 (Pl-Sarajedini). These 
observations were taken as a part of a large survey of 
65 globular clusters, and consist of 4 deep exposures and 
one shallow exposure in each of F606W and F814W. The 
observations were taken on 22 Jun 2006, so the baseline 
between this set and the first is 4.07 years. The right 
panel of Figure [1] shows that this field is centered on the 
cluster, and on the GO-9442 mosaic. 

In the following section, we describe how we con- 
structed an extensive catalog based on the i?F435W and 
-Rf625W exposures of the 10' x 10' GO-9442 data set; 
then we cross-identify the same stars in the other data 
sets and measure proper motions. 

3. REDUCTIONS 

In this section we describe how we constructed a ref- 
erence frame and a catalog of photometry and proper 
motions. In the subsequent sections, we will then use 
this catalog to determine the cluster center and study 
the surface-density and velocity-dispersion profiles. 

3.1. Reference frame 

The first task in the analysis is to construct a reference 
frame. We began by measuring all the bright, isolated, 
unsaturated stars in all the deep exposures of GO-9442. 
There are 36 exposures total in H-alpha, and 27 in each 
of -Bf435W and i?F425W- To do this, we used the publicly 
available software program img2xym_WFC . 09x10 (Ander- 
son & King 2006), which employs an empirically con- 
structed, spatially variable library PSF to determine a 
position and a flux for each star in each exposure. The 
positions are corrected for distortion using the correc- 
tions in Anderson (2005). 

We began by analyzing the nine major pointings inde- 



3 

pendently. For each, we took the central i?F625W observa- 
tion as the basis for the reference frame of that pointing, 
linearly transformed the star positions from all the other 
images of that pointing into that frame, and determined 
an average position for each star that was found in at 
least three exposures. 

The next step was to stitch together the nine overlap- 
ping frames. There was only minimal overlap between 
the frames. We did not want to allow for linear terms 
in the frame-to-frame transformation, since small errors 
in such terms could lead to large errors in the extrapo- 
lated field, so we stitched the pointings together by solv- 
ing only for offsets and rotations, assuming the frames 
to have the same scales and off-axis linear terms. We 
first added the directly adjacent side fields to the central 
frame, than added the corner fields to this plus-shaped 
intermediate frame. 

To improve this reference frame, we found a linear 
transformation from each exposure into the new frame, 
based on the positions of common stars. For each star 
in the reference frame, we thus had between 7 and 40 
estimates for its position (depending on how many im- 
ages overlapped at that point), and we averaged these 
positions together to improve the reference frame. After 
a few such iterations, we found the remaining stitching 
errors to be less than 0.01 pixel. 

The above reference frame was constructed in a 
distortion-corrected frame based on the central pointing. 
To align our frame with RA and Dec, we found the lin- 
ear transformation from this frame into the frame con- 
structed in Anderson ct al. (2008; hereafter AOS), based 
on the GO-10775 Treasury data, which has a pixel scale 
of 50 mas/pixel and was aligned with RA and Dec, and 
has the Harris (1996) cluster center at [3000,3000]. (We 
will refer to this as the "traditional" center, since it is 
based on the literature at the time.) We conformally 
transformed our reference frame into this frame and 
added 4000 to each coordinate, so that our entire field 
would be contained within the range [1:14000,1:14000], 
with the center close to [7000,7000]. We refer to this 
frame as the "master" frame. 

At this point we made a digression, to put our reference 
frame on an absolute basis, even though it is not needed 
in the present study. The guide-star catalog, and hence 
the header information, can contain astrometric errors 
of up to 1.5 arcsecond (Koekemoer et al. 2005), so we 
cross-identified over 6000 of our stars with the 2MASS 
catalog in order to construct a more accurate absolute 
astrometric calibration of our reference frame. We found 
that the 2MASS frame was offset by 0.5 arcsecond in 
declination and 2 arcseconds in right ascension from the 
frame given in the GO-10775 drz-frame header. We also 
found that there was a small orientation difference (0.1°) 
between the drz-based frame and the 2MASS frame. The 
WCS (world-coordinate system) parameters in the fits 
headers of the stacked images that we provide were con- 
structed to match as well as possible the positions of stars 
in our frame with those in the 2MASS catalog. Skrut- 
skie et al. (2006) argue that the 2MASS positions for 
well-measured stars should be good to about 50 mas in a 
random sense, and good to 15 mas in a systematic sense. 

We emphasize that our proper motions do not depend 
greatly on the quality of our astrometric solution. Each 
star will be measured relative to its nearby neighbors us- 
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ing only individual flat-fielded (-fit) images; the refer- 
ence frame contributes only to our scale and orientation. 
Nonetheless, for future use of the catalogs provided here 
it is still desirable to have the best possible reference 
frame. 

3.2. Image stack 

The transformations from the frame of each exposure 
into the master frame allow us to generate stacked im- 
ages of the field in this frame. Our images were generated 
in a manner akin to drizzle with a pixf rac parameter 
of (see Fruchter & Hook 2002). These stacked images 
provide a simple representation of the field that allows 
us to independently evaluate how our star-finding algo- 
rithms performed (our finding was done on the individual 
images, not on the stacks). 

The stacked images are 14000x14000 pixels, covering 
11.6 arcminutes on a side, and were constructed for the 
F435W, F625W, and F658N filters of GO-9442. We pro- 
vide these images as a part of the data release with this 
paper. 

3.3. Constructing the Catalog 

Once we had determined a reference frame, we could 
construct a star catalog in this frame. To do this, we 
used the same software program that was used for the 
reductions for the 65 clusters in GO-10775. The details 
of this program can be found in AOS. Briefly, the routine 
reads in the available images (shallow and deep) for two 
different filters (-Bf435W and i?F625W, in this case) and 
auxiliary information that allows it to map each exposure 
into the reference frame, both astrometrically and pho- 
tometrically. It then goes through the reference frame in 
patches that are 25x25 pixels in size, and identifies all 
the stars in each patch, measures them, and records their 
parameters in a file. The routine is able to deal well with 
short and deep exposures; if a star is saturated in the 
deep exposures, it is found in the short exposures where 
possible. It also finds stars iteratively, first finding the 
bright stars, then removing them to find the fainter stars 
that might not stand out distinctly in the im-subtracted 
images. Finally, the routine is robust against identifying 
PSF artifacts or diffraction spikes as stars. 

The GO-9442 data set differs from the typical GO- 
10775 cluster data set. Here, the master frame was 
14000x14000 WFC pixels in size, and consisted of a 
mosaic of pointings, whereas in the GO-10775 set, the 
observations consisted of medium-sized dithers about 
one central pointing and could fit comfortably within a 
6000 X 6000-pixel region. The typical point in the field 
under study here had coverage of 3 B and 3 R observa- 
tions. In order to be included in the catalog, we insisted 
that a star be found in at least 4 out of these 6 ob- 
servations independently; if there were fewer exposures 
available at a given point in the field, the criteria were 
relaxed accordingly (see AOS). 

The automated finding program identified 1,164,317 
stars in the i3F435W and -Rf625W images. We plotted the 
star lists on top of the stacked images to verify that the 
routine had identified all the stars we expected it to. This 
inspection confirmed that there were very few stars that 
were missed by the automated routine. The few missed 
stars fall into two classes: medium-brightness stars in the 
vicinity of saturated stars, and blended pairs of stars that 




Fig. 2. — A stack of the deep F625W images at the center 
of the cluster. The yellow circles identify the stars found in the 
S435 +-R625 list. The red-circled stars were found in the -ffess list, 
but not the -B435 + -R625 list. The large green circle is the location 
we have identified for the center (see Section |4]|. The field shown 
is 21" X 16". 

had one star that was brighter in i?F435W and the other 
brighter in i?F625W- The former were passed over by 
the program because they were too likely to be artifacts 
near saturated stars, and it was better to err on the side 
of caution. The latter were missed because we required 
each star to be found as a dominant peak in at least four 
independent images. Both of these loss mechanisms can 
be quantified with artificial-star tests (see below). 

To get an idea of how many real stars were missed, 
we ran the same finding program on the F65SN images 
(which had less saturation and four exposures at each 
pointing), requiring a star to be found in at least 3 out 
of 4 exposures, and cross-identified the stars in the two 
catalogs. There were about 30,000 stars (3% of the to- 
tal) that were found in the H-alpha images that were 
not found in the -Bf435W and i?F625W observations. We 
looked at where these missing objects were located, and 
they were almost exclusively found in the situations de- 
tailed above. 

In Figure [H we show the central region of the clus- 
ter. The stars found jointly in i3F435W and i?F625W are 
marked in yellow. The stars that were found in i?F658N 
(H-alpha) but not in -Bf435W and -Rf435W are marked in 
red. Of course it would be good to have a list that con- 
tains every single star in the cluster, but a compromise 
must always be made between including the marginal 
star and including image artifacts as stars. There will 
always be stars that are missed; the important thing is 
to be able to quantify this incompleteness. 

3.4. Calibrating the photometry 

The fluxes for the stars measured in the above routine 
are reported in instrumental magnitudes, — 2.51ogio^, 
where Z is the scaling that matches the effective PSF 
model to the stellar profile. The PSF model was normal- 
ized to have a flux of 1.00 within a radius of ten _f It 
pixels. The short exposures have been zero-pointed to 
match up with the deep exposures, therefore the task 
of calibration simply involves determining the zero-point 
that will bring the instrumental photometry of the deep 
exposures into the VEGAMAG system. 
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Fig. 3. — The CMD of all the stars found in the automated 
search of the F435W and F625W exposures of GO-9442. The level 
where saturation sets in for the deep exposures is indicated, as 
well as the number of stars that are saturated. Most of the stars 
above the deep-image saturation line were measured well in the 
short exposures. Saturation in the short exposures sets in at about 
4 magnitudes brighter than this. These particularly bright stars 
were measured by fitting the PSF to the unsaturated part of the 
star's profile. Note the multiple main-sequence turn-offs, which 
have been well-discussed in the literature (i.e., Ferraro et al. 2004 
and Villanova et al. 2007). 

Since the GO-9442 data set has already been cah- 
brated by Villanova et al. (2007), we simply adopted 
the photometric zero points from that project. They are 
2'Pf435W = 32.043, ^Pf625W = 31.991, and ZPpessN = 
28.904. Figm-e [3] shows our (-Bf435W — ^F625w) color- 
magnitude diagram (CMD) in both the instrumental and 
in the calibrated systems. We will continue to report 
many of the analyses here in the instrumental system, 
since it is easier to assess errors in terms of signal-to- 
noise ratio when working directly in photo-electrons. For 
reference, saturation in the deep images sets in at an in- 
strumental magnitude of around —13.75, and a star with 
a magnitude of —10.00 has a signal-to-noise ratio of 100. 
In the rest of the paper, when referring to the instru- 
mental system, we will report magnitudes as mF435W: 
and when referring to the calibrated VEGAMAG pho- 
tometry, we will refer to -Bf435W- 

3.5. Artificial- star tests 

The GO-9442 images are not terribly crowded. There 
are about 1.2 million stars in 14Kxl4K pixels, thus the 
typical separation is about 10 pixels. Nevertheless, the 
fact that there are quite a few extremely bright stars 
makes it impossible to find all the faint stars. And even 
though the density changes by less than a factor of four 
from the center to the corners of the field, it is still im- 
portant to have an idea of the completeness for each 
brightness level of star at different places within the field. 
Therefore, we decided to run artificial-star tests in order 
to gauge the incompleteness and measurement quality. 

The mechanism for artificial-star (AS) tests is de- 




FlG. 4. — Central swath of the stacked image (15" by 15" in each 
panel). On the left, the yellow circles indicate stars that were found 
between —10.5 and —11.5 (in F435W instrumental magnitudes). 
On the right we show where stars in the same brightness range 
were successfully recovered from our artificial-star tests (green), 
and where they were not recovered (red). The center (determined 
in Section |4]l is shown in both panels with a blue l"-radius circle 
around it. 

scribed thoroughly in AOS. Briefly, since our finding soft- 
ware operates on one small patch of the field at a time, 
we can afford to do artificial-star tests in serial, one at a 
time, rather than doing them in parallel in many batches. 
One great benefit of this is that artificial stars never in- 
terfere with each other, so by throwing many stars in 
successively we can simulate throwing in a high density 
of them. This makes it easy to do a detailed study of the 
completeness in the vicinity of the center. 

We performed 500,000 individual artificial-star tests. 
For each, we chose a random F435W instrumental mag- 
nitude between —14 (saturation) and —5, and chose the 
F625W magnitude that placed the star along the fiducial 
main sequence. The artificial stars were inserted with a 
flat distribution in radius, so that we would probe the 
central regions more than the outer regions. For exam- 
ple, we ended up inserting over 700 stars within the cen- 
tral arcsecond in radius, and over 51,000 stars within the 
central 10 arcseconds. 

As a demonstration of how the AS tests worked, in 
Figure [4] we show the same central field as in Figure [2l 
and the locations where the artificial stars were and were 
not recovered. In this region, we observed 261 real stars 
between —10.5 and —11.5 in F435W instrumental magni- 
tudes. We ran 3847 artificial-star tests in this magnitude 
range over this 225-square-arcsecond region, and recov- 
ered 3402 of them, resulting in a completeness of 88 %. 

Figure [5] shows the results of the artificial-star tests for 
different radial bins as a function of instrumental mag- 
nitude. It is clear that the completeness at the center 
is well over 90% for the brightest stars, and is over 75% 
throughout the cluster for stars with instrumental mag- 
nitudes brighter than —10, which corresponds to stars 
about 4 magnitudes below the upper SGB. 

3.6. Measuring the Proper Motions 

Figure [1] shows that proper motions can be measured 
in two fields: one field that is largely centered on the 
cluster and another that is east south-east of the center 
by about 4 arcminutes. Table [T] lists the observations we 
have available for each field. 

Our procedure for measuring proper motions was to 
measure the stars individually in each exposure for each 
epoch, and then to combine the many independent mea- 
surements to construct proper motions for the stars that 
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Fig. 5. — The completeness as a function of magnitude for stars 
at different distances from the cluster center. 

can be found in multiple epochs. The agreement among 
the many independent observations for each star at each 
epoch gives us a handle on the errors in the positions, 
and hence in the displacements and proper motions. 

The GO-10252 images are well matched depth-wise to 
the GO-9442 images, thus we should expect the same 
stars to be measured equally well in both epochs. By 
contrast, the GO-10775 images, which cover the central 
field, are a factor of four shallower than the GO-9442 data 
set, and therefore have different signal-to- noise ratios for 
the same stars. 

3.6.1. Measuring the star positions 

The first step in determining proper motions is to mea- 
sure an accurate position for each star in each individual 
exposure at each epoch. To do this, we once again ran 
the star-measuring program img2xym_WFC . 09x10 to con- 
struct a list of sources in each exposure, using empirically 
determined PSFs with a spatially-constant perturbation 
to account for breathing-related focus changes from ex- 
posure to exposure. 

The routine produced a list of positions and instru- 
mental magnitudes for the reasonably bright and isolated 
stars identified in each exposure. Since our aim was to in- 
clude only the unsaturated stars that could be measured 
accurately, we adopted finding parameters such that we 
would find unsaturated stars with at least 250 counts in 
their brightest four pixels and with no brighter neighbors 
within 4 pixels. We also confined our attention to the 
deep exposures for each filter/epoch combination, since 
there is generally only one short exposure at each loca- 
tion and it is hard to get a good handle on astrometric 
errors without multiple exposures. Finally, we treated 
the central and adjacent major-axis proper-motion fields 
separately, since there is only marginal overlap between 
them. 

3.6.2. The central field 



For the central field, we have 10 exposures in the first 
epoch in filters -Bf435W, -Rf625W5 and -ffpessw, and 8 ex- 
posures in the second epoch in filters Vpeoew and /f814W • 
We reduced each exposure as described above, produc- 
ing a list of X, y, and m (instrumental magnitude) for 
each source. We treated the two WFC chips for each 
exposure independently, so we had a total of 36 lists of 
stars, one for each chip in each exposure. We corrected 
the positions for distortion and cross-identified each list 
with the reference-frame list, allowing us to develop a 
least-squares linear transformation from each exposure 
into the reference frame. Since the central field covers 
only a fraction of the GO-9442 field, we restricted our- 
selves here the subset of 422,561 stars contained within 
the region [4000:9999,4000:9999]. 

The above linear transformations yielded an estimate 
for the reference-frame position for each star from each 
exposure in which it was found. These global transforma- 
tions can contain some small systematic residuals, due to 
variations in the distortion solution caused by breathing 
or other effects. The consequence of these errors is that 
the transformed positions for the stars in one part of the 
reference frame are all shifted by some (small) amount. 
In order to remove these residual transformation artifacts 
from the data, we determined a local adjustment for each 
observation of each star as follows. We first identified all 
the neighboring stars within a 100-pixel radius that were 
unsaturated and had a S/N > 40 (excluding the target 
star itself). For each exposure where the star was ob- 
served, we then found a robust average offset between 
the globally transformed positions of the neighbor stars 
for that exposure and the average reference-frame po- 
sitions for the stars. This average offset provided the 
correction to the transformed position of the target star 
for that exposure. 

Figure [6] shows the residuals before and after this cor- 
rection for the first image for each epoch. The typical ad- 
justment is 0.01 pixel, which is about the accuracy that 
we expect for the static distortion solution (Anderson 
2003). But not all of this adjustment comes from static 
distortion-solution error; much of it clearly varies from 
exposure to exposure due to breathing-related changes in 
focus. The effect of this correction (and, indeed, of the 
star-based transformations themselves) is that the mo- 
tion of each star is measured relative to that of the av- 
erage motion of its neighbors. Since outliers are rejected 
in the process of determining the average motion of the 
neighbors, this means that each star's motion is com- 
puted relative to the bulk cluster motion at its location. 
Section 13.6.41 discusses what this calibration procedure 
means in practice for the motions that we measure. 

At this point, we had in hand corrected reference- frame 
positions for a large number of stars measured in a large 
number of independent exposures. We combined the 
data in several ways. First, for each filter/epoch com- 
bination, we averaged all the data together to arrive at 
an average position, magnitude, and RMS error for each 
star. We also combined the astrometric data together 
by epoch. (Since the F658N data have a lower signal-to- 
noise ratio at the faint end, we include F658N astrome- 
try only when the star had an instrumental magnitude 
of i?F658N < -10, indicating a S/N of 100). 

The proper-motion results for the central field are il- 
lustrated in Figure [71 Each data set consisted of 3 to 
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Fig. 6. — (Left) In the top panel, we show the x-residual between 
the first F435W image in GO-9442 (j61p05vuq) and the reference 
frame for a 100-pixel-tall strip through the cluster center. Since 
the reference frame is based on images taken at that epoch, these 
residuals should show just distortion errors and should have no 
proper motions. The middle panel shows the local "boresight" 
correction used for each star. The bottom panel shows the residuals 
after the local adjustment. (Right) Same, but for the first exposure 
in the second epoch (j919a7ftq). These residuals contain both 
distortion errors and actual proper motions, and therefore show 
more scatter than in the left panels. 
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7. — Panels (a) through (c) show the distribution of 

T^, the total two-dimensional error in the position for the 

first-epoch F435W, F625W, and F658N data, respectively. This 
error was computed by taking the RMS of the scatter among the 
multiple observations of the master-frame position for each star 
and then dividing by the square root of the number of observa- 
tions. The dotted line drawn at 0.02 pixel indicates the boundary 
between the well and poorly measured stars. Panels (d) and (e) 
show, respectively, the same for the F606W and F814W data of 
the second epoch. 



4 observations of each star. We constructed an average 
reference-frame position for each star in each data set, 
then computed the error in this position from the RMS 
about this average. The two-dimensional error in the av- 
erage position at a single epoch is typically 0.01 pixel. 
This error was achieved for the data set for each filter at 
each epoch. 

The proper motions were constructed by subtracting 
the first-epoch positions from the second-epoch positions 
and dividing by the time baseline. The multiple data sets 
we had at each epoch allowed us to construct two inde- 
pendent observations of the proper motions. We con- 
structed one displacement comparing the second-epoch 
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Fig. 8. — Panel (a) compares two independent measurements of 
the X displacement between epoch-2 and epoch-1. The agreement 
between the two is an indication of the quality of the proper mo- 
tions. Panel (b) shows the same for the y displacement Panel (c) 
shows the resulting 2-d displacement distribution, after combining 
the data from all filters for each epoch. Panel (d) zooms out to 
show the cluster and field stars (points beyond 0.25 pixel are dis- 
played as crosses, for clarity). The time baseline between these 
epochs is 4.07 years. The pixel scale is 50 mas/pixel, and x-axis 
of the reference frame is aligned with the negative right-ascension 
axis. 

F606W positions with the first-epoch F435W positions, 
and a second displacement by comparing the second- 
epoch F814W with the first-epoch F625W. The top pan- 
els of Figure [8] compare these two independent measure- 
ments. The distribution along the 45°-direction repre- 
sents the actual motion measured, found to be the same 
in both halves. The distribution in the 135°-direction 
represents the error in our measurements. The typical 
displacement between epochs is about 0.15 pixel. The 
displacement of the typical star is measured with a frac- 
tional error of less than 10% (from the aspect ratio of 
distribution in the upper panels). 

We combined all the data for each epoch and con- 
structed a single proper motion for each star, which we 
plot in the lower panels as vector-point diagrams. There 
are very few outlier stars, meaning that we have done a 
good job selecting stars with small errors. It is worth 
pointing out that our finding criteria have not selected 
against high-motion stars: the field stars in this diagram 
demonstrate that we are finding stars with displacements 
of up to 1 pixel. Thus, we are clearly sensitive to clus- 
ter members with more than 10 times the typical motion. 
The baseline for this central field is 4.07 years, thus a dis- 
placement of 0.08 pixel corresponds to about 0.02 pixel 
per year in each coordinate, or 1 mas/yr. 

We measured proper motions for the 108,507 stars that 
could be found in at least two images of each epoch, with 
an RMS of position measurements within each epoch of 
less than 0.03 pixel (per coordinate). In addition, we 
defined a "high-quality subset" of the data that should 
have more uniform PM errors by selecting those stars for 
which this position RMS was less than 0.02 pixel and for 
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Fig. 9. — Left: all the stars from the catalog in the region of 
overlap between the central pointing of GO-9442 and GO-10775. 
Right: the stars found to have reliable proper motions. On the 
right of each CMD, we add up the number of stars in each 1- 
magnitude bin. 

which the instrumental m.F435w photometry was brighter 
than —11. Figure [H] shows a CMD of the stars that qual- 
ify for the high-quality subset. Our high-quality motions 
are limited to stars on or below the sub-giant branch 
(SGB) and those brighter than tof435W = ^H, a few 
magnitudes below the turnoff. Stars much brighter than 
the SGB were saturated in the deep GO-10775 exposures. 
Most of the proper-motion analysis in this paper and Pa- 
per II will use this high-quality sample of 53,382 stars. 
Within this sample, the typical proper-motion error is 
better than 0.1 mas/yr in each coordinate, corresponding 
to 0.008 pixel over the 4-year baseline. Over the magni- 
tude range where we could measure reliable motions, we 
measured between 80% and 90% of all the stars. 

3.6.3. The adjacent major-axis field 

In a similar way to that described above, we con- 
structed proper motions for the adjacent field along the 
major axis. Since the GO-10252 data straddled two of 
the GO-9442 pointings, we had to use two first-epoch 
pointings here, which doubled the number of exposures 
to handle, but the overall reduction was generally the 
same. To qualify for this proper-motion catalog, we re- 
quired that a star be found in at least 4 first- and 4 
second-epoch exposures. (Because of possible distortion 
errors at the edges of the out-of-focus GO-10252 data set, 
we insisted on a star being found in at least 4 exposures 
to ensure that the internal errors will be able to indicate 
when stars have distortion issues.) 

We found 61,293 stars for which we could measure 
proper motions in the adjacent field. As above for the 
central field, we defined a high-quality sample of stars 
that should have uniformly good motions. To qualify for 
this sample, a star had to have an RMS of measure- 
ments within each epoch of less than 0.02 pixel, and 
had to be brighter than r7iF435w — —11- The high- 
quality proper-motion catalog for the adjacent field con- 



tains 19,593 stars. 

3.6.4. Proper Motion Zero-Points 

The proper motions we have measured here are not ab- 
solute proper motions, as is generally true for other glob- 
ular clusters observed with HST (see McLaughlin et al. 
2006 for a broader discussion). HST can measure abso- 
lute proper motions only when it has absolute reference 
points to measure displacements against. The Omega 
Gen fields are so crowded that there are no detectable 
galaxies that could serve as absolute reference points. 
Also, the stability and repeatability of HST (while su- 
perb) are not sufficient to compare different frames di- 
rectly. For example, there are small changes in focus 
due to temperature-induced telescope breathing. These 
cause scale changes that need to be calibrated out. More- 
over, the nominal orientation (roll) repeatability of HST 
is 0.003 degrees. If left uncalibrated, this would induce 
apparent solid-body rotation of the cluster that exceeds 
the internal motions. Hence, to be able to derive relative 
proper motions, we need to use general six-parameter 
linear transformations to match exposures to each other. 

The application of linear transformations implies that 
scale, rotation, two skew terms, and two translation 
terms cannot be measured from the data. The removal 
of scale and skew terms is not an issue, because the dy- 
namics of a cluster in equilibrium do not display time- 
variability in these terms. The translation terms imply 
that we can't measure the mean motion of Omega Cen 
with respect to the Sun. This does not affect the inter- 
nal cluster dynamics, which is our primary concern here. 
Therefore, the only important fact to keep in mind is that 
all proper motions reported here are measured modulo an 
undetermined solid-body rotation component. We show 
in Section 5.2.3 of Paper II that this does not impact our 
modeling. Ground-based line-of-sight studies of Omega 
Gen reported by van de Ven et al. (2006) show that any 
true solid-body rotation in our central field should be 
quite negligible, even though we can't measure it di- 
rectly here. Note in this context that any unaccounted 
for solid-body rotation component would generally tend 
to decrease the IMBH mass implied by models, since it 
would act to lower the kinematical gradient between the 
inner and outer parts of the cluster. 

To be able to remove small time-variations in high- 
order geometric distortion terms, we have gone one step 
further in our calibration. Namely, we measured each 
star relative to its local neighbors (see Section I3.6.2p . 
All the neighbors used should be moving with the clus- 
ter, but each will have some random internal motion. We 
measured each star against hundreds of its neighbors, so 
the random internal motions should average out. There- 
fore, the proper-motion zeropoint is different for each 
star, and equals the mean projected motion of cluster 
stars at its position on the sky. In other words, all mean 
motion is removed, not just solid-body rotation. 

The "local-neighbor correction" procedure allows the 
most accurate measurement of the internal proper- 
motion dispersions for the cluster. However, this cor- 
rection is applied by choice, and not by necessity. We 
would not want it to remove an important signal that is 
in fact present in the data. Mean motion in the radial 
direction should not exist in an equilibrium cluster. How- 
ever, there could be differential rotation (i.e., the part of 
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Fig. 10. — This figure exhibits the average local corrections 
applied when deriving the proper motions, as function of position. 
The average vectors were derived as described in the text for 18 X 
18 bins that cover the inner 1'.5 X 1'.5 . Vectors that extend to the 
edge of the circle correspond to a displacement difference of 0.01 
pixel, which is typical of geometric distortion residuals. The size of 
the corrections is too small (approximately one eighth of the RMS 
cluster motion) to have any dynamical influence. The morphology 
of the vector field indicates that we are not inadvertently removing 
rotation intrinsic to Omega Cen. 



the rotation curve that is not a hnear function of radius). 
A Keplerian rise in rotation around an IMBH would fall 
in this category. There is no evidence for such differen- 
tial rotation in any existing data, but we can verify this 
directly with our HST data. Without the local-neighbor 
corrections, the data are fully sensitive to small-scale in- 
ternal rotations and motions. 

Figure [10] visualizes the size of the local-neighbor cor- 
rections that were applied. Our set of stars with good 
proper motions were distilled into boxes 100 pixels on 
a side. We computed the proper motion for each star 
two ways. First, we took the esimate generated above, 
which involved removing from each position measured 
in each image the average of the local neighbors. Sec- 
ond, we skipped this "boresight" correction, using only 
6-parameter global linear transformations to relate posi- 
tions in each image to the master frame. The difference 
between these two measured positions is shown as the 
vector in the plot. A vector that reaches the edge of the 
circle corresponds to a displacement difference of 0.01 
pixel. 

The vectors in Figure [TOl correspond to the mean local 
motions that were subtracted. Two things are evident. 
First, the corrections are very small. The typical ad- 
justment is less than 0.01 pixel, which is characteristic 
for geometric-distortion residuals. This corresponds to 
about an eighth of the observed proper-motion disper- 
sion, and is too small to have dynamical importance for 
the cluster. Second, there is no overall rotation pattern 
visible. The correction field is patchy and incoherent. 
Again, this is as expected for geometric-distortion resid- 
uals, and this cannot be due to internal motions of the 



cluster. More specifically, there is no excess rotation near 
the cluster center that would be typical for Keplerian ro- 
tation around an IMBH. 

As a further check, we also used the the methods of 
Paper H to determine the rotation curve (mean tangen- 
tial proper motion as function of radius) for the proper- 
motion catalogs with and without the local correction. 
The difference in rotation curves was less than 1 km/s at 
all radii in the central field. 

In summary, our catalog formally contains proper mo- 
tions after subtraction of all mean motion in the plane 
of the sky. However, detailed analysis and ground-based 
data indicate that any mean motion that the cluster may 
in fact possess is dynamically irrelevant at the central- 
field radii of our HST study. 

In Sections [6] and U\ we present an initial analysis of 
the proper motions. Paper II will combine these motions 
with existing proper-motion, radial- velocity, and surface- 
profile data from the literature and will construct com- 
prehensive dynamical models for the cluster. 

3.7. Catalog and Image products 

In the previous sections we described how the WFC im- 
ages were analyzed to construct star and proper-motion 
catalogs. These catalogs are available along with the on- 
line version of this paper. 

Specifically, we have analyzed the 10' x 10' field of 
dataset GO-9442, which goes out to beyond 2 core radii. 
We have one file (DATA . CATALOG . XYVI) that contains 
1,164,328 stars found from running our automated find- 
ing algorithm on the F435W and F675W images^. Ta- 
ble [2] gives the column- by-column description of this file. 
These measurements can be converted to calibrated pho- 
tometry using the zero points included in the table (from 
Section [3^ . We also provide the fits images referred 
to in Section[221 which contain a WCS header and which 
are in the same coordinate system as the star catalogs. 

We also ran 500,000 artificial-star tests and subjected 
them to the same finding and measuring algorithm as 
was used on the real stars. These tests are reported in 
file DATA . ART . XYVI, which is also available in its entirety 
with the on-line article. The artificial-star tests report 
the same measurement quantities as in Table [5] for the 
real stars; but for the AS tests, we also have the input po- 
sitions as well, recorded in DATA. ART. INPUT. It is worth 
pointing out that just because something was "found" in 
an artificial-star test, that does not mean that the input 
star was actually recovered; it is important to compare 
the input and output positions and fluxes to determine 
that the inserted star, and not a brighter pre-existing 
neighbor, was indeed recovered. 

Finally, we provide the proper-motion measurements 
in separate files for the central and major-axis fields, 
DATA . PMs . CEN and DATA . PMs . MA J. The central field has 
motions for 108,507 stars and the adjacent major-axis 
field for 61,293 stars. Table [3] gives a column-by-column 
description of the proper-motion data files. The "good" 

^ Upon careful inspection of our bright stars, we found 112 
bright stars spread throughout the field in the initial catalog that 
were actually detector artifacts along the bleed-pixels of extremely 
bright saturated stars. We identified these stars easily as they 
had no counterpart at all in the H-alpha or the short images, thus 
we have removed them from our catalog. This is the only non- 
automated aspect of our catalog. 
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TABLE 2 

The columns in the star catalog, DATA. CATALOG. XYVI. 

Columns 3 and 4 give instrumental magnitudes, which can 

BE transformed TO THE VEGAMAG SCALE BY ADDITION OF 
THE LISTED ZEROPOINTS ZPvEGA- STARS BRIGHTER THAN ABOUT 
— 13.75 ARE SATURATED IN THE DEEP EXPOSURES, AND THEREFORE 
HAVE LARGER SYSTEMATIC ERRORS. We REPORT THE RMS 
SCATTER IN INDIVIDUAL OBSERVATIONS OF EACH STAR. In 

principle, the errors in columns 1 through 4 can be 
constructed by dividing these rmss by the square root of 
the number of observations (columns 10, 12 or their sum). 

However, not all errors are random, so this may 
underestimate the errors. the rms can be used to help 
select stars with good measurement consistency. 



Column 


Description 


1 


X position in the reference frame 


2 


y position in the reference frame 


3 


'"■F435W1 instrumental B magnitude, ZPvega = 32.043 


4 


"iF625Wi instrumental R magnitude, ZPyEGA = 31.991 


5 


RMS scatter of the single-exposure x positions 


6 


RMS scatter of the single-exposure y positions 


7 


RMS scatter of the single-exposure r(iF435W observations 


8 


RMS scatter of the single-exposure m,F625W observations 


9 


number of B images where star could have been found 


10 


number of B images where star was found well 


11 


number of R images where star could have been found 


12 


number of R images where star was found well 



TABLE 3 

The COLUMNS in the proper-motion CATALOGS, DATA.PMs.CEN 

AND DATA.PMs.MAJ. The PROPER-MOTION ERRORS ARE 
DETERMINED BY ADDING THE ERRORS IN THE AVERAGE FOR EACH 

EPOCH IN QUADRATURE, ASSUMING NO INTER-EPOCH 
TRANSFORMATION ERROR. THE MAJOR/mINOR AXIS PROJECTIONS 

were made assuming the major axis to be 100° east of 
North. 



Column 


Description 


1 


X position in the reference frame 


2 


y position in the reference frame 


3 


Bf435W, calibrated VEGAMAG B magnitude 


4 


i?F625W 1 calibrated VEGAMAG R magnitude 


5 


^x, proper motion in x, mas/yr 


6 


^ly, proper motion in y, mas/yr 


7 


0"mua,, proper-motion error in x, mas/yr 


8 


'^muy, proper-motion error in y, mas/yr 


9 


A'maji proper motion along major axis 


10 


/^mini proper motion along minor axis 


11 


A^el , number of exposures with good positions in epoch 1 


12 


Ng2} number of exposures with good positions in epoch 2 


13 


g, the "good" flag, indicating small internal errors 



flag refers to stars that are part of the high-quahty sub- 
set defined in Section I3.6.2[ which is used in the later 
sections of this paper and in Paper II. 

In the following sections we will use these catalogs to 
determine the center of the cluster (§ H]), to extract a 
surface-density profile (§ [5]), and to search for a cusp 
in the velocity distribution that could be indicative of an 
IMBH at the center (§[5]). Finally, in §[71 we will compare 
the proper motions for stars in different populations. 

4. DETERMINING THE CENTER 



In general, the center of a cluster can be determined 
to a precision of about a/^/N, where a is roughly the 
core radius and describes the fall-off of the spatial distri- 
bution, and N is the number of stars used in the center 
determination. For a cluster like Omega Cen, which has 
a core radius of 2'.5 (150"), this means we would need 
over 150^ = 25, 000 stars to measure a center to within 
an arcsecond. These stars would have to be spread out 
to well beyond a core radius to give us a good handle 
on the center. (The center is defined as much by where 
stars are not as by where they are.) 

In the case of a Gaussian distribution, the most accu- 
rate way to measure the center is in fact to take a simple 
centroid of the star positions. However, there are sev- 
eral factors that can complicate this simple solution. For 
example, if the region surveyed is not symmetric with re- 
spect to the center, then we must include some correction 
for this effect on the centroid. Similarly, if the list of stars 
used in the determination is not complete throughout 
the field, then we must consider how this will affect the 
center determination. The presence of bright giants can 
prevent us from finding the more plentiful main-sequence 
stars in their vicinity, and the small-number statistics of 
the giants' distribution can result in patchy incomplete- 
ness. All these limitations can be overcome; but it is 
important to address the issues carefully. 

In this section we employ several different ways of de- 
termining the center from our data. We determine the 
density center from star counts, and the kinematic center 
from proper motions. In both cases we use two separate 
methods to find the point of symmetry: one based on 
contours and the other based on pie-slices. As an in- 
dependent cross-validation, we also derive the center of 
unresolved cluster light in a 2MASS mosaic image of the 
cluster. We find that all methods yield centers that agree 
to within the calculated errors. We compare these centers 
with previous determinations available in the literature. 

4.1. Iso-Density Contour Centroids 

The first way in which we determined the center was by 
constructing iso-density contours, and fitting ellipses to 
them. We began by sifting the stars in our 14000 x 14000- 
pixel reference frame into an array of 140x140 bins, 
where each bin is 100x100 pixels (5" x 5"). We found 
that these bins did not contain enough stars to generate 
smooth contours, so we subsequently ovcrbinned the dis- 
tribution, making each bin correspond to 500 x 500 pixels, 
or (25" X 25"). 

The upper-left panel of Figure [TT] shows the contours 
for the stars with tof435w < — 9. The edges of the field 
are clear, as is the region north of (9000, 8000) where 
there is a gap in the coverage (see Figure [T]) Our aim 
was to fit ellipses to the valid parts of these contours. The 
contours traced in black are the ones that we trusted to 
be far enough from the edges of the field to be valid. 

In panel (B) we show the points that fall along every 
third contour (for clarity), and the ellipse that we fit- 
ted to each. In these fits, along each contour we used 
only the points that were well away from the gap region. 
Panel (C) shows the ellipses fit to all the contours. The 
parameters of the fitted ellipses are given in Table |4l 

Panel (D) in Figure [TT] shows the centers of the con- 
tours shown in Panel (C). The average position for the 
center is (6727,6810) with an RMS of about 45 pixels in 
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Fig. 11. — Panel A: The isodensity contours for all stars 
with mF435M' < ~9- Contours are separated by 100 stars (per 
25" X 25" box). The blue contours are contaminated by the edges 
of the field and are not used. Panel B: The points along contours 
700, 1000, 1300, 1600, 1900, 2200, and 2500, and the ellipses that 
were fit to them. Points in the gap-region in the northwest were 
not used. Panel C: The fitted ellipses (the parameters are given in 
Table|4]l. Panel D: Distribution of the ellipse centers. The circle is 
centered on the average and has a radius equal to the RMS of ob- 
servations about the average. This RMS is also listed in the panel. 
The uncertainity in the average position is smaller by \/7V, where 
JV > 4 is the number of independent contours. 



TABLE 4 

Summary of the fitted-ellipse parameters for the 

CONTOURS from STARS BRIGHTER THAN ^435 = —9. THE DENSITY 

is reported in stars per 25" x 25" bin; the fitted center is 
pixels in the master frame: a is the semi-major axis and b 
the semi-minor axis, in wfc pixels; p. a. is the position 
angle of the major- axis, measured in degrees e from n. 
The last column represents the ratio of the minor to the 

MAJOR axis. 



Density 


Fitted center 


H^ + A) 


P.A. 


B/A 


700 


6780 6905 


5055 


100 


0.87 


800 


6750 6810 


4648 


108 


0.86 


900 


6695 6810 


4320 


104 


0.86 


1000 


6740 6785 


4018 


102 


0.86 


1100 


6720 6825 


3726 


106 


0.87 


1200 


6740 6790 


3471 


106 


0.88 


1300 


6760 6755 


3207 


104 


0.88 


1400 


6765 6785 


2982 


102 


0.89 


1500 


6780 6805 


2749 


102 


0.89 


1600 


6775 6800 


2539 


96 


0.89 


1700 


6720 6840 


2303 


92 


0.90 


1900 


6755 6815 


1841 


102 


0.87 


2000 


6670 6795 


1618 


118 


0.80 


2100 


6625 6780 


1387 


118 


0.86 


2200 


6620 6845 


1082 


134 


0.83 


2300 


6685 6835 


784 


132 


0.94 


2400 


6730 6840 


526 


56 


0.87 


2500 


6760 6740 


262 


102 


0.54 



each coordinate. Since the contours are not entirely in- 



dependent of each other, we cannot simply construct an 
error in the center by dividing the RMS by the square- 
root of the number of contours fit. However, since we 
have at least four independent contours (based on their 
separation and the size of our bins), the center should be 
good to about 20 pixels, or 1". 

4.2. The Pie-Slice Method 

The focus of the previous approach was to find the 
center of the cluster by identifying the centers of vari- 
ous iso-density contours of the stellar distribution. In 
this section, we describe how we divided the cluster into 
pie slices and determined the center by finding the loca- 
tion about which the stars are most symmetrically dis- 
tributed. 

The complication for this second approach is that the 
broad PSF halos of the few bright giants can cause a 
patchy incompleteness in the distribution of the faint 
stars, which can in turn cause us to find more stars on one 
side of the cluster than on another, and thus to misiden- 
tify the center. The contour-based approach in the pre- 
vious section skirted this issue somewhat by focusing on 
the outer contours, where crowding is lower and the gra- 
dient is steeper; but the pie-slice approach will need to 
deal with incompleteness all the way into the center. 

In order to measure an unbiased center, we had to con- 
struct a star list that is free of any incompleteness that 
would bias the center. This list must also contain enough 
stars to give us a statistically reliable result, so we could 
not simply use the few bright stars, which are essen- 
tially 100% complete. We clearly need to include the 
more plentiful faint stars, but including them requires 
taking careful account of incompleteness. Our artificial- 
star tests from Section 13.51 provided one way to assess 
the incompleteness in the field. Unfortunately, it was 
not feasible to run enough artificial-star tests to assess 
completeness on scales small enough to adequately sam- 
ple the local region around each bright star. What we 
needed was a simple map of the field that could tell us 
how bright a star must be in order to be found in a 
particular place in the field. If we had such a map, we 
could then create a mask that is symmetric about any 
presumed center, so that the exclusion zones caused by 
the bright stars would be symmetric, and hence would 
not bias the center determination. 

4.2.1. Completeness Mask 

Our goal in constructing a mask was to determine, for 
a particular bright star, where in its vicinity a given faint 
star could be reliably found. To answer this question, we 
identified the bright stars with mF435w < — 16 (about 
10 X saturation) and looked at the distribution of their 
found neighbors as a function of distance from them. 
This distribution is shown in Figure [121 At a distance 
of 10 pixels, only stars that are within 5 magnitudes of 
the bright star's fiux can be found; fainter stars than this 
cannot be found reliably. Out to a distance of 20 pixels, 
stars up to 7.5 magnitudes fainter than the bright star 
can be reliably found. This d versus Am relationship is 
quantified by the red line drawn in the figure. A few 
stars above this line are found, but essentially, all the 
stars below this line are found. 

Note that there are very few artifacts above this line; 
normally bumps in the PSF at regular distances and off- 
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Fig. 12. — The distribution of found stars in the vicinity of 
brighter stars {mp^^^iY < —16). The sohd line (red in the on- 
Hne version) indicates the faintest magnitude of star that could be 
definitively found at a given distance from a brighter star. 

sets are identified as stars by automated finding algo- 
rithms. We have almost no such artifacts in our lists 
because our initial finding procedure had a model of the 
PSF out to tens of pixels and the procedure was careful 
not to identify faint stars where appreciable PSF features 
are likely to be present (see Fig. 3 of AOS). If we had 
identified every bump in every image as a star, then there 
would be a large number of detections in a distinct pat- 
tern above the line. (Fig. 1 of McLaughlin et al. 2006 
shows what the un-screened PSF artifacts look like in 
this space.) 

The red line in Figure [Til allows us to construct a zone 
around each bright star that tells us which stars stand 
out clearly at each distance. We brought the information 
from all the bright stars across the entire field together 
in order to construct an image that tells us how bright 
a star must be at each point in the field in order to be 
definitively found. The value of each pixel in this image 
tells us the magnitude of the faintest star that could be 
reliably found at each point. Figure [T3l shows the F435W 
image of the central region of the cluster. The middle 
panel of the figure shows the finding mask constructed 
from these stars. 

This mask essentially corresponds to an additional se- 
lection requirement: stars must have been identified by 
our automated procedure and they must also satisfy this 
location-based requirement. The goal of the mask is not 
to perfectly reflect our artificial-star tests, but rather to 
provide an estimate of which stars could definitively be 
found at various points in the field. We validated this as- 
sumption using the artificial-star tests: we verified that 
artificial stars were indeed recovered essentially every- 
where the mask says they should be recovered. 

One final note about the finding mask. There is a gap 
in the north-west part of our field (see Fig[T]). Since no 
stars could be found where there was no coverage, we 




Fig. 13.— On the left we show a 10" X 10" F435W image of 
the of central field. In the middle panel, we show the mask for the 
same region. The mask shown here has been gray-scaled to show 
the areas where stars of brightness —8, —10, and —12 can and 
cannot be found. The mF435w ^ —8 stars can be reliably found 
only in the white areas. The mF435W ^ —12 stars can be found 
everywhere except for the very darkest areas. The actual mask 
used for the analysis has much more gradation than is shown here. 
On the right, we show the mask after it has been symmetrized 
about a presumed center (very close to the actual center). 

assigned the mask a value of —25 in this region. For ref- 
erence, the brightest star in the field has an instrumental 
magnitude of -20.28 in F625W. 

4.2.2. Symmetrization 

The inferred mask made it possible to create star lists 
that have symmetric incompleteness properties with re- 
spect to any adopted center. With this tool in hand, we 
next determined the best center using a pie-slice proce- 
dure similar to that used in McLaughlin et al. (2006). To 
do this, we assumed an array of trial centers. About each 
trial center we divided the stars into pie slices centered 
upon that location and computed a statistic to compare 
the cumulative radial distributions of the stars in oppos- 
ing pie slices. We finally identified the center as the place 
in the field that has the most similar radial distributions 
in the opposing pie slices. 

Specifically, we decided to use eight pie slices, arrayed 
in cardinal and semi-cardinal directions (see Figure [H)). 
We explored centers within the range of reference-frame 
coordinates [6400:7200,6400:7200], with a trial center ev- 
ery 20 pixels (1") in each coordinate. In order to en- 
sure that our star lists did not have any asymmetric bi- 
ases, for each prospective center (/c, Jc), we generated 
a list of stars by asking whether each star in the cata- 
log would be found both at its own location in the mask 
(i* and at the location on the other side of the center 
(/c — [i*— /c], Jc — [j* — Jc\)- If the star satisfied both cri- 
teria, it qualified for the center-determination list. The 
rightmost panel in Figure[T3]shows what the mask looked 
like when symmetrized about a point near the ultimate 
center. A star that could not be found at one place in 
one slice would be excluded from the corresponding place 
in the opposing slice. 

We next selected the stars that both satisfied this 
symmetrized-mask criterion and were brighter than 
"^F435W = —9 within 4000 pixels radial distance (200"), 
and distributed them among the eight pie slices accord- 
ing to their azimuthal angle. We sorted the stars in 
each slice by distance from the adopted center, then for 
each of the four opposing pie-slice pairs we computed 
E„=mm(jVi,jV2) V^i"^) - '"2(")l> which is the integrated 
difference between the radial distance r versus cumula- 
tive A'cnci distribution functions for the opposing slices. 
This statistic was constructed for the four opposing pairs 
of slices at each trial center. 
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Fig. 14. — The two contour plots to the left show the goodncss- 
of-center parameter for the pie slices shown below. The symbol at 
the center represents the paraboloid fit to the 9x9 points centered 
on each minimum. In the upper right, we compare the innermost 
contours and the fitted centers. In the lower right, we add the two 
goodness-of-fit metrics to get an overall best fit. The filled circle is 
the center from NGB08, and the 6-pointed star is the center from 
the contour-based determination. The small square and triangle 
correspond to the center of the cardinal and semi-cardinal contours, 
respectively. 



TABLE 5 

Summary of center determinations using the pie-slice 
METHOD. All positions are reported in our reference 

FRAME. 



Selection Number Cardinal Semi-cardinal Together 

B435 < -9 235,000 (6718,6820) (6718,6820) (6703,6801) 

S435 < -10 205,000 (6723,6837) (6738,6784) (6723,6811) 

B435 < -11 150,000 (6736,6854) (6759,6807) (6747,6829) 



'71F435W < — 11, and found very similar results. Table [5] 
summarizes our findings. Overall, the centers determined 
by the different brightness cutoffs are consistent. 

The center we found from the contour-based study 
in Section 14.11 was (6727,6810), which is in excellent 
agreement with the (6723,6811) location found here for 
mF435w < —10. In what follows we will adopt a center of 



(6725,6810) for the number density of Omega Cen, with 
an estimated error of ~1" in each coordinate. 

4.3. Kinematic Center from the Proper Motion 
Dispersion Field 

An independent method to verify the results from the 
previous subsections is to find the kinematic center. In 
the present context, that means identification of the 
point of symmetry in the proper-motion velocity disper- 
sion field. Unlike the star-count analysis, the kinematic 
center determination has the advantage of being inde- 
pendent of any incompleteness corrections (so long as 
the kinematics of a star do not affect whether or not a 
star makes it into the catalog). 

We took the full proper-motion catalog for the central 
field and binned the stars onto a grid of 1" x 1" pixels. 
The few stars believed to be outside the cluster (as in- 
dicated by Figure [8]i) were excluded from the analysis. 
For each pixel we calculated the number of stars as well 
as the second proper- motion moment (/x^ -I- fiy) summed 
over the stars. This yields two images, which we will 
call Nq{x, y) and N2{x, y). These images are quite noisy, 
due to the relatively small number of stars per pixel. We 
therefore applied a two-dimensional top-hat smoothing 
kernel to each image to increase the signal-to-noise ratio. 
This yields images N{){x,y) and N2{x,y). The implied 
smoothed image of the one-dimensional RMS proper- 
motion a{x,y) is given by [N2/{2No)Y/^{x,y). In prac- 
tice we found that a top-hat kernel radius of 25" yields 
adequate results. 

The resulting RMS proper-motion image (j{x, y) is 
shown in Figure 1151 There is a well-defined symmetric 
distribution around a broad central peak. The irregular 
outer boundary of the image represents the intersection 
of the two (approximately square) ACS/WFC fields for 
epochs 1 and 2, respectively. The cluster is not centered 
within this boundary, due to the particular details of 
the observational pointings. The closest boundary line 
is ~ 70" South from the cluster center. In the kernel 
smoothing we excluded all pixels that are outside the 
boundary region of the catalog. Nonetheless, properties 
of the map within 25" from the boundary may be some- 
what affected by artifacts induced by the absence of data 
outside the boundary. To ensure that this would not bias 
the determination of the cluster center, we restricted our 
analysis of the map to radii R < 45" from the cluster 
center. 

We performed ellipse fits to the image to determine 
the symmetry point of the map. We adopted the average 
ellipse center for semi-major axis lengths between 35"- 
45" as our final estimate. The result lies at (A^^, Aj,) = 
(—1.1", —2.0") from the adopted star-count center. To 
determine the random errors in this center determina- 
tion we performed Monte-Carlo tests. We created pseudo 
data sets by populating the same area of the sky covered 
by our catalog, with the same number of stars. The stars 
were drawn from the projected number density profile de- 
rived in Paper II. Each star was assigned proper motions 
in the x- and y-directions by drawing random Gaussian 
deviates from distributions of dispersion a. The value 
of cr was chosen to be a slowly decreasing function of 
radius, consistent with our measurements in Section [5] 
and Paper II. Proper motion errors were added based on 



Figure [U] shows the result for our pie-slice analysis. 
The four pairs of opposing slices allow us to construct 
two independent estimates of the center, shown in the 
contours in the top plots. We fit a paraboloid to the 
central 9x9 points in the contour plot to arrive at the 
best-fit center for each panel. The plot in the upper-right 
panel shows the agreement between the two centers, and 
the two central contours. The plot on the bottom-right 
shows the contours constructed from all 4 pairs of op- 
posing slices. The difference between the two center es- 
timates is 2", indicating that our average center is prob- 
ably good to about 1". In this determination of the cen- 
ter, we used about 235,000 stars. With a core radius of 
150", we would expect the center to be accurate to about 
150" /a/235000 or about 0.3 arcsecond; this is not a huge 
discrepancy if we consider that the cluster distribution is 
not Gaussian and the star lists used here did not extend 
out beyond two core radii. 

We conducted the same procedure with a subset of our 
stars, using only those stars with mF435w < — 10 and 
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Fig. 15. — Velocity field of the one-dimensional proper-motion 
dispersion, smoothed with a 25" r adiu s top-hat kernel. The orien- 
tation of the image is as in Figure 1111 The irregular outer border 
of the image is the intersection of the two (approximately square) 
ACS/WFC fields for epochs 1 and 2, respectively. At the chosen 
level of smoothing, the gradient in dispersion between the center 
and the edge of the field corresponds to ~ 2 km/s. Ellipses are 
fits to the contours of the map, with semi-major axis lengths of 
30"-70". The lower dot (blue) is the center of the ellipse (shown 
in blue) that has a semi- major axis length of 40". This is used as 
an estimate of the HST kinematic center, i.e., the symmetry point 
of the map. The small circle is the 68% confidence contour around 
the kinematic center, as determined from Monte-Carlo simulations. 
The confidence contour shows that the kinematic center is consis- 
tent with the HST star count center (upper dot, black). The center 
adopted by NGB08 (black square) is 12" away, and is inconsistent 
with both of the HST centers. 

Gaussian deviates from randomly selected error bars in 
the observed catalog. Each pseudo data set was analyzed 
in similar fashion as the real catalog. The Monte-Carlo 
kinematic centers thus determined had an RMS scatter 
of 2.9" per coordinate around the input center, with no 
bias. The elongation of the a{x, y) contours in Figure [T5l 
was found not to be statistically significant. EUipticities 
as large as the observed value (at R — 40") of ~ 0.23 
happened by chance in 14% of the Monte-Carlo simula- 
tions. 

To further ensure the robustness of the results, we also 
experimented with alternative approaches. These used 
percentiles of the proper-motion distribution (instead of 
the RMS), different grid sizes, different smoothing ker- 
nels, or different kernel sizes. The results were always 
consistent with those quoted above. We also applied 
the method to the high-quality proper motion subsam- 
ple, instead of the full proper motion catalog. This too 
yielded consistent results, albeit with larger uncertain- 
ties.^ This is because the full catalog provides a larger 
sample size than the high-quality subset. Although the 
full catalog has somewhat larger proper motion uncer- 
tainties, these uncertainties have essentially no impact 
on the analysis (although they are fully included in the 

^ Specific results were included in an earlier preprint version of 
this paper. These are now superseded by the analyses presented 
here. 



simulations). They are always much smaller than the 
cluster dispersion. Therefore, the errors in the map of 
Figure [15] are determined primarily by the number of 
stars that contribute to each pixel (Act w a /\/2N), and 
not by the individual per-star proper motion uncertain- 
ties. For the adopted kernel size, N k, 6000 near the 
center. This yields random errors of ~ 0.15 km/s, which 
is much smaller than the gradient in the dispersion map 
(in fact, the size of the kernel was purposely chosen to 
make this the case). 

Our method is purposely designed to measure the sym- 
metry point of the proper motion map. One could use 
the peak of the map, but this quantity is much more af- 
fected by shot noise, due to the intrinsically low spatial 
gradients in the core of Omega Cen. Specifically, the 
peak of the proper motion dispersion map in Figure [15] 
is at (Ax, Ay) — (—6", —7") from the adopted star-count 
center. However, the Monte-Carlo simulations show that 
the peak pixel, when used as an estimate of the cluster 
center, has an RMS uncertainty of 6.7" per coordinate. 
This uncertainty gets even larger if one chooses a smaller 
smoothing kernel size than the 25" used here. So while 
the observed peak is statistically consistent with the sym- 
metry point, it is a much more unreliable estimate of the 
true cluster center. 

This makes an important point that underlies all our 
analyses: Omega Cen has such a large and (nearly) ho- 
mogeneous core, that to determine its true center most 
accurately, one must locate the symmetry point of its 
large-scale distribution. Measures of possible small-scale 
peaks in density, star light, or kinematics, are much more 
susceptible to shot noise. While this does not make such 
estimates incorrect, they can only be interpreted if their 
uncertainties are rigorously quantified. 

4.4. Kinematic Center from the Pie-Slice Method 

As an alternative approach to determination of the 
kinematic center, we also used a pie-slice method. For 
a given trial trial center, we adopt a polar grid with 
N azimuthal wedges {N chosen to be even) and M ra- 
dial bins along each wedge. The radial bins are linearly 
spaced and have size S. The value of M is chosen so 
that M X S 100", and the grid therefore encompasses 
most of our central-field catalog. For each bin we use our 
full proper motion catalog to calculate the proper motion 
dispersion ct and its error Act as in Section [4.31 For each 
pair of radially opposing bins we calculate the quantity 
(cti — CT2)^/(ActJ -I- Act|), which measures the extent to 
which the dispersions in the two bins are statistically 
consistent with each other. When summed over all the 
opposing bins, this yields a quantity. The number of 
degrees of freedom TVdf is normally M x N /2. However, 
we exclude pairs that do not have a sufficient number of 
stars in each bin to yield a meaningful dispersion. Each 
pair thus excluded reduces A'df by one. We map the 
quantity A^^ = — A'^df on a grid of trial centers with 
1" spacing. We smooth the resulting map with a Gaus- 
sian with a 2" dispersion (this smoothing is not required, 
but was found to reduce the uncertainties in the resulting 
kinematic center estimates by ~ 15%). The trial center 
with the lowest Ax^ is the position with respect to which 
the proper motion dispersion field is most radially sym- 
metric. 

To minimize any dependence of the results on the 
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somewhat arbitrary choices of N and M, we repeated 
this procedure for N = 4,6,8,10,12,14,16 and S = 
6, 9, 12, 15, 18, 21, 24 arcsec. We smiimed the Ax^ maps 
from all the 49 combination of N and S to obtain one 
grand-total Ax^ map. The lowest Ax^ in this map oc- 
curred at (Aa;,Aj,) = (-0.4", -1.5") from the adopted 
star-count center. We used the same pseudo-data sets 
as in Section 14.31 to characterize the uncertainties in the 
result. Each pscudo data set was analyzed in similar 
fashion as the real catalog. The Monte-Carlo results had 
an RMS scatter of 1.9" in x and 2.4" in y. This mea- 
sures the uncertainties in the inferred kinematic center. 
The Monte-Carlo results were found to be slightly bi- 
ased, and were centered on {Ax,Ay) = (—0.4", —0.6"). 
Since this bias can be calculated and corrected, it is 
does not affect the accuracy of the final result. The 
bias, and also the fact that the scatters differ in the 
two coordinate directions, are due to the specific geom- 
etry of our proper motion dataset.'^ Upon bias correc- 
tion, the pie-slice estimate for the kinematic center is 
{A^, Ay) = (0.0" ± 1.9", -0.9" ± 2.4"). We repeated the 
whole analysis also with only the high-quality observed 
proper motion sample. This again yielded consistent re- 
sults, but with somewhat larger uncertainties. 

The kinematic center thus determined using pie-slices 
is very consistent with the value derived from contour 
fits in Section 14.31 In fact, all methods that we have 
explored to determine the kinematic center yielded the 
same answer to within the uncertainties. The pie-slice 
method yields the lowest uncertainties, because it uses 
more of the large-radius information in the catalog (dat- 
apoints with 45" < R < 100" are now included). We 
therefore use the pie-slice result as our final estimate for 
the kinematic center. 

Our analysis provides the first time that the kinematic 
center of any globular cluster has been accurately deter- 
mined. The kinematic center of Omega Cen was found to 
be consistent with the star-count center. This is what is 
expected in an equilibrium system, and therefore should 
not come as a surprise. However, this is important since 
it provides an independent verification of the center po- 
sition determined from the star-count analysis. 

4.5. Comparison with the NGC08 Center and Other 
Literature Values 

In ^3.11 we used the 2MASS point-source catalog to 
tie our master reference frame to RA and Dec. The cen- 
ter we identified in this frame corresponds to (a, 5) = 
(13:26:47.24, -47:28:46.45), with an error of about an 
arcsecond in each direction. 

We were unable to find any other centers for Omega 
Cen in the literature that had quoted errors. The most 
recent determination of the center was done by NGB08, 
who found the center to be at (13:26:46.04,-47:28.44.8). 
It turns out that this absolute position was measured 
in the drizzled ACS images, and contains some system- 
atic error due to errors in the guide-star-catalog positions 
(Koekemoer et al. 2005). If we identify their center from 
the star field shown in their Figure 2 (which corresponds 

^ The bias can be avoided by using the method only with data 
closer to the center, instead of going out to 100" where the az- 
imuthal coverage becomes incomplete. However, that would pro- 
duce larger uncertainties. 
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Cen. Estimates of the cluster center from Harris (1996; plus), van 
Leeuwen et al. (2000; diamond), and NGB08 (cross) are shown 
in blue. None ol^ these authors provided error estimates for their 
centers. The green boxes identify the two fields studied with IFU 
spectroscopy by NGB08, one of which was believed to include the 
cluster center. The circles/ellipses mark the cluster center positions 
determined in the present paper using vari ous i ndependent meth- 
ods, namely HST star counts (red; Section 14.21 1. HST proper mo- 
tions ( mage nta; Section |4.4| I. and 2MASS unresolved light (cyan; 
Section 14.61 1. The sizes of the circles/ellipses indicate the 68.3% 
confidence regions of the estimates (note: for a two-dimensional 
circular Gaussian probability distribution, the circle that encloses 
68.3% of the probability has a radius that equals 1.516 times the 
one-dimensional error bar). Our estimates are mutually consistent, 
and their uncertainties clearly rule out the prev ious ly reported val- 
ues in the literature (see discussion in Sections 14. 71 and 15.31 1. 

closely to the absolute coordinate they report from the 
WCS header of image j61p05weq_drz), then their cen- 
ter corresponds to (6962,6881) in our master frame, or 
(13:26:46.08, -47:28:42.9) using the 2MASS astromet- 
ric reference frame, which is (-1-237,-1-71) pixels or 12.3 
arcseconds away from our center. 

FigurelUlprovides in finding-chart format the locations 
of our star-count and kinematic centers determined from 
our HST analysis, with the error estimates indicated, as 
well as the center and IFU fields used by NGB08. NGB08 
do not estimate the error in their center, so it is hard to 
say that our centers are formally in disagreement. How- 
ever, given that they felt their center was within their 
5" X 5" IFU field, it seems safe to assume our centers are 
in significant disagreement. We explore this further in 
Section [5?3l 

Harris (1996) reports (13:26.45.9, -47:28:37) for the 
center in his on-line catalog, van Leeuwen et al. (2000) 
found a center of (13:26:45.756, -47:28:42.780) using the 
positions of the giant stars brighter than i? = 16. This 
center is adopted by van de Ven et al. (2006) and Castel- 
lani et al. (2008). The Harris and van Leeuwen centers 
are also indicated in Figure [HI They are closer to the 
NGB08 center than to the center positions inferred by us. 
We explore the differences between these ground-based 
centers and our HST center further in Section [4.71 
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Fig. 17. — Upper left: 2MASS J-band mosaic centered on Omega 
Cen. Upper right: Result of filtering the image to derive the un- 
derlying background (as described in the text). Lower left: Result 
after additional top-hat smoothing, with contours drawn in. Lower 
right: Close-up of the yellow-boxed region in the upper-left panel. 
Locations of the contour centers are shown in cyan. The red circle 
is the HST star-count center derived here and adopted in our mod- 
eling; the 1" radius of the circle corresponds to the one-dimensional 
error bar. The NGB08 center is shown in blue. 

In Table [6l we summarize the centers from the htera- 
ture and compare them with the center determined here. 

4.6. Center from 2MASS 

Because of the disagreement between the center po- 
sitions derived by us and those previously published 
in the literature, we sought non-HST-based data that 
could serve as an independent cross-check. The pub- 
licly available data from 2MASS proved suitable for this 
purpose. We downloaded a wide 2MASS image mosaic 
of the cluster using the NVO/ISRA mosaic service at 
[littp://hachi.ipac.caltech.cdu:8080/montage/ . The up- 
per left panel of Figure [17] shows the J-band mosaic im- 
age. The pixel scale is one arcsecond per pixel, and the 
field covers 14' x 11'. 

Since the 2MASS image was taken at near-IR wave- 
lengths, it is extremely sensitive to the red-giant stars. 
In order to minimize the shot noise from bright stars on 
our analysis, we applied a "filter" to the image to isolate 
the underlying unresolved light. The filtering operation 
consisted of going through the mosaic image pixel by 
pixel, examining the surrounding pixels within a radius 
of 20 pixels (20"), and computing the tenth percentile of 
this neighbor distribution. The image on the top right of 
Figure flTl shows the result. It shows the general profile 
of the background light, but it is still patchy on account 
of the few bright stars. 

We next smoothed this background-filtered image with 
a circular top- hat kernel that had a radius of 20 pixels, 
and arrived at the image in the lower left panel of Fig- 
ure [iTl This image is smoother and allows the derivation 
of a contour-based center. The blue lines are image con- 
tours between 2 core-radii and 0.75 core radii, drawn at 
equal intervals of intensity. The adopted radial range is 
well-suited for a determination of the center. It avoids 
the very central region, where the spatial gradient is too 
low to yield high accuracy. It also avoids, the outer re- 
gion, where the results can be biased by uncertainties in 
the 2MASS sky-background subtraction process. 

We fitted ellipses to the contours as in ii4.11 and show 
the resulting centers overplotted in cyan in the bottom 
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right panel of Figure [iTl This panel shows a close-up of 
the yellow-boxed region in the upper-left panel. The av- 
erage of the centers lies at (Aj,, Ay) = (2.1", 1.9") from 
the adopted HST star-count center. The RMS scatter 
in the ellipse centers is 2.1" per coordinate. Given the 
kernel size for the percentile-filter and blurring opera- 
tions, we estimate that about a third of the 12 contours 
between 2 core-radii and 0.75 core radii are statistically 
independent. Therefore, the error in the mean position 
per coordinate is ^ j^fi = 1.0". Given that the er- 
ror in the HST center is also ^ 1", we find the centers 
determined from HST and 2MASS to be in acceptable 
statistical agreement. By contrast, the 2MASS position 
is 9.8" from the NGB08 center, and even further from 
the Harris (1996) and van Leeuwen et al. (2000) centers. 
It is therefore inconsistent with those centers. 

Use of unresolved light, especially in the near-IR, is 
more prone to possible systematic errors than our analy- 
ses based on HST star counts and kinematics. Therefore, 
we would not assign the 2MASS result the same level of 
confidence as our HST results, despite the similar ran- 
dom error. However, the 2MASS result is useful as a 
cross-check. The fact that it agrees with our HST results 
indicates there is no reason to suspect some fundamental 
problem with using either HST or ground-based data to 
determine the cluster center. 

4.7. Understanding the differences in ground-based 
centers 

To understand how so many previous investigations 
could infer centers that differ so much from the cen- 
ters found here, we need to address two separate ques- 
tions: (1) why did previous authors who used ground- 
based data, such as Harris (1996) and van Leeuwen et 
al. (2000), derive centers that differ so significantly from 
that derived here from 2MASS?; and (2) why did NGB08, 
who used some of the same HST data that we have, de- 
rive a center that differs so significantly from our star- 
count and kinematic-center results? We address the first 
question here. The second question is closely tied to the 
determination of cluster number density profile from the 
HST data, and we therefore discuss it in Section 15.31 be- 
low. 

We simulated a ground-based image by combining the 
R-band photometry from our (nearly complete) HST cat- 
alog with a broad, 3-arcsecond FWHM PSF. This is 
shown in the upper left panel of Figure [TH We also 
generated an image that represents the number counts 
for stars with S/N > 100 in the deep HST F625W expo- 
sures, distilled into the same 0.5"-pixels as the simulated 
image. These medium-brightness stars should not suf- 
fer much from incompleteness. The star-count image is 
shown in the upper right panel of Figure 1181 

It is clear in the upper left panel that there are lines 
of bright stars in both dimensions that generate bright- 
ness enhancements that cross at the rough location of 
the "traditional center". The latter is at the center of 
the box, roughly 15" W and 10" N of the HST center de- 
rived above, which is marked by the red circle. It is also 
clear that the visual centroid of the bright-star distribu- 
tion is N and W of the HST center that we have identi- 
fied. At the same time, the upper right panel shows that 
the bright, S/N > 100 sources are much more evenly dis- 
tributed. They appear centered, as expected, around the 
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Fig. 18. — Upper left: simulated image with 3" FWHM res- 
olution of Omega Cen using our R catalog; Upper right: number 
density of S/N > 100 stars from our catalog; Lower left: horizontal 
profile of the two data sets; Lower right: vertical profile of the two 
data sets. The arrows delimit the regions used in the creation of 
the profiles. The red circle is the HST star-count center derived 
here and adopted in our modeling; the 1" radius of the circle corre- 
sponds to the one-dimensional error bar. The "traditional center" 
(Harris 1996) adopted in many studies of Omega Cen is at the 
center of the inner box. The offset shown on the horizontal axis of 
the bottom panels is measured with respect to the HST star-count 
center. 



previously derived HST center. Even by eye it is clearly 
evident that the S/N > 100 sources are not centered on 
the center of the box. The dark areas of the distribution 
extend more into the left and bottom parts of the panel 
than they do into the right and top parts. 

We extracted horizontal and vertical profiles from each 
of these images. The lower left plot shows the horizontal 
profile across the 30"-tall swaths (between the horizontal 
arrows). The lower right plot shows the same profile for 
the vertical direction. The blue points correspond to the 
simulated R-band image, and the green points to the 
number counts. It is clear that the traditional center 
corresponds to a location where there is a coincidental 
concentration of bright stars along both the x and y axes. 

It is beyond the scope of this paper to determine how 
each previous determination of the center went wrong. In 
fact, many of the previous determinations did not pro- 
vide a catalog or method description of sufficient detail to 
allow easy reexamination. Also, none of them provided 
an error estimate. Nevertheless, we have shown that a 
few bright stars are responsible for a light enhancement 
that is not coincident with the overall density peak of the 
more numerous, but fainter, stars. This light enhance- 
ment is close to the positions where previous authors 
estimated the cluster center to be. So it is likely that 
the methods previously used on ground-based data were 
disproportionately influenced by this light enhancement, 
and that this led to the historical mis-identification of 
the center. 

The analysis in Figure [TSl does not imply that it is im- 



possible to derive the cluster center from ground-based 
data. After all, we did manage to derive an accurate 
center from 2MASS data. However, what is clear from 
these analyses is that at a minimum one must adopt spe- 
cial methods to mitigate the shot noise from bright stars 
(e.g., through filtering as used in our 2MASS analysis). 

5. SURFACE-DENSITY PROFILE 

Our nearly complete catalog of stars, covering the in- 
ner two core radii of the cluster, allows us to calculate a 
definitive radial density profile for the inner part of the 
cluster. Note that the entire cluster cannot be described 
by a single radial profile, since the cluster contains stars 
of different masses at different stages of relaxation. If 
we consider only the profile of the evolved stars, then we 
will suffer small-number statistics, since only about 1% 
of the stars in our catalog are above the SGB. Our strat- 
egy here has been to measure the number-density profile 
for stars within a range of magnitudes. We note that 
number-density profiles do a much better job describing 
the star distribution than do surface-brightness profiles. 
A star at the tip of the RGB is 100 times as bright as 
a star on the SGB, yet the two stars have essentially 
the same mass and should both be equally good trac- 
ers of the density distribution. Since the brighter stars 
are not better tracers of the star density, it makes sense 
that we should not give them more weight. Furthermore, 
since there are many times fewer evolved stars than SGB 
and main-sequence stars, surface-brightness profiles suf- 
fer much more from shot noise than do number-density 
profiles. 

5.1. Profile from the new HST Catalog 

Our strategy in computing the number-density profile 
was straightforward. We divided the field into concen- 
tric annuli centered on the center derived in Section |4l 
Within each annulus, we determined the number of stars 
observed in each magnitude bin. We also determined 
the completeness for the magnitude bin for that annulus, 
based on the artificial-star tests. We then constructed 
an average surface density within each annulus for each 
magnitude bin by dividing the number of found stars by 
the completeness. 

Figure [19] shows the surface-density profile for three 
different brightness ranges, from the stars brighter than 
?7iF435W = — 13 (just above the turnoff), to stars just be- 
low the turnoff at mF435w to stars well down the 
main sequence at 7tif435w ^ ~9. It is hard to know 
exactly which profile we are most interested in. The 
profiles for the brightest stars tell us about the most 
massive luminous stars, but the fainter-star profiles have 
more stars, and therefore provide better constraints on 
the profile. All the profiles appear reasonably well fit by 
the single-mass King model we provide as a reference. 
Paper II fits the data with parameterized models, which 
will serve as the basis for our dynamical modeling. 

There is no evidence of a sharp rise at the center in any 
of the profiles. Indeed, the central bin, which extends out 
to about 2" in radius, has fewer stars than the surround- 
ing bins. Section 2] showed that our center is accurate to 
about 1", so the center should be well contained within 
our central bin. The completeness for all these bins is 
greater than 80%, even at the very center. Looking at 
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Fig. 19. — In each panel, we show the surface-density profile for 
stars in the labeled magnitude range (in stars per square arcsccond, 
with a linear scale). For reference, the bottom of the RGB is 
roughly at mp^^e^yv = —13.5. The filled circles have been corrected 
for completeness. The open circles (green in the on-line edition) 
correspond to the density of stars actually observed (for clarity, we 
plot this only for the inner 50"). The error bars indicate the V-^obs 
errors. The radial bins are 2 arcseconds wide; for the innermost 5 
bins we note at the top N^hsy the actual number of stars observed. 
The curve (red in the on-line version) corresponds to a single-mass 
King model with a core radius of 2'.5 and a tidal radius of 59', 
shifted vertically to fit each profile. A more detailed quantitative 
analysis of model fits to the number density profile is presented in 
Paper II. 

the image in Figure[2l it is clear that there is not a signif- 
icant increase in the vicinity of the center. Quantitative 
constraints on the central slope of the number density 
profile will be presented in Paper II. 

5.2. Comparison with the NGB08 Surface-Brightness 

Profile 

In their paper, NGB08 measured the snviace- brightness 
profile, rather than the suvi&ce- density profile. They 
found that it rose with a power law with a logarithmic 
slope of —0.08 ± 0.03, which they deemed to be in sig- 
nificant disagreement with a flat core. This is at odds 
with our result, in that we found a different center and 
a relatively flat profile about it. 

Since our data set is a super-set of the data they used, 
we decided to try to reconstruct their profile using our 
star catalog. In their work, NGB08 measured bi-weights 
of the background light distribution in the actual F435W 
image they used. We approximated their calculation by 
listing all the pixels within each of their annuli (centered 
on their center) and taking the 25th percentile of the 
points. This is plotted as the filled black points in Panel 
(A) of Figure [20l and it agrees quite closely with the 
profile shown in Figure 1 of NGB08. 

The NGB08 aim in measuring the surface brightness 
from the background was to be less sensitive to the bright 
giants and more sensitive to the numerous fainter stars. 
Unfortunately, most of the light in the background comes 
not from a large population of unresolved stars, but 



ACTUAL IMAGE 
SIMULATED IMAGE 
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Fig. 20. — In Panel (a) we show that if we use the NGB08 radial 
bins, and a percentile-based estimate of the background, we observe 
the same power-law radial trend they saw (filled black points). 
The open points represent the same procedure, but performed on 
a simulated image, as described in the text, which contains no faint 
unresolved stars. Panel (b) shows the number-density profile about 
their center. Panel (c) shows the catalog stars about the NGB08 
center. 



rather from the PSF halos of the bright stars. This can 
be seen from a simple examination of the luminosity func- 
tion. In an external galaxy where we cannot detect stars 
below the SGB, the luminosity function (LF) is seen to 
increase steeply from the brightest stars to well below 
the detection limit, and as such there is a large reservoir 
of faint stars just below the detection limit. By contrast, 
in a globular cluster such as Omega Cen, we resolve and 
count a large fraction of all the stars that are present. 
The stars that are not individually resolved contribute 
almost no light. In this sense, the analysis of unresolved 
light from HST data is different than for ground-based 
data. For example, the 2MASS images analyzed in Sec- 
tion 14.61 do not resolve individual stars on or below the 
SGB. So unhkc the HST case, in the 2MASS images there 
is a large repository of unresolved stars that contribute 
a significant amount of light. 

As a consequence of the fact that the LF is relatively 
flat below the detection limit, most of the background in 
the HST images comes from the bright-star halos rather 
than faint, undetected stars. Measuring the profile of the 
HST image background is therefore equivalent to measur- 
ing the profile from a blurred version of the bright stars. 
To demonstrate this, we simulated the field by taking the 
positions and brightnesses of all the stars in our catalog 
and using PSFs that go out to 100 pixels (5"), based on 
the encircled-energy curves in Sirianni et al. (2005). We 
performed the same 25th-percentile-based procedure on 
our simulated image, and arrived at the open dots, which 
trace the actual observations extremely well. There ap- 
pears to be a constant offset of ~100 electrons between 
the real and the simulated profiles, which can easily be 
accounted for by diffuse light scattered by even more than 
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100 pixels — the bright core of this cluster extends out to 
3000 pixels (2.5 arcminutes). 

In Panel (b) we show the number- density profile about 
the center adopted by NGB08. This profile does not 
increase as steeply and monotonically as the surface- 
brightness profile. To understand qualitatively how the 
surface-brightness profile was found to be more cuspy 
at the NGB08 center, Panel (c) shows the placement 
of the NGB08 apertures in the field. Their central 50- 
pixel annulus happened to just contain a large number 
of turnoff-brightness stars. However, there is no evident 
central concentration of the more plentiful fainter stars, 
as a true cusp would suggest. The second annulus (from 
50 to 100 pixels) contains a large number of bright RGB 
stars. The halos of these encircling stars may help to 
explain why the background in the inner annulus was 
observed to be elevated. 

5.3. Understanding the differences in HST-based 
centers 

Both our study and the NGB08 study used HST data 
to determine the center of Omega Gen. It is therefore of 
interest to examine in more detail the differences between 
our results. 

In Section 0] we have used the distribution of three 
quantities to determine the center of Omega Gen, namely 
HST star counts, HST proper motions, and unresolved 
2MASS background light. For each quantity we used a 
contour-based method, and for the HST-derived quanti- 
ties we also used a pie-slice method. The five different 
analyses all give the same answer to within the uncertain- 
ties of 1-2 arcsec. This effectively rules out the possibility 
that unexplored bias in any of the methods might have 
significantly affected the analysis. What the methods 
all have in common is that they identify the symmetry 
point of the cluster using data that extends significantly 
from the cluster center. This naturally uses the full size 
of the data set to reduce shot noise. The center thus 
identified is 12.3" away from the NGB08 center. Hence, 
the NGB08 center is ruled out as the symmetry point 
of the cluster at 5-lQa confidence by each of 3 different 
quantities, analyzed with several different methods. 

NGB08 used a pie-slice approach on HST star count 
data to estimate the cluster center. The general method 
they used is described in Noyola & Gebhardt (2006), but 
few details are provided about the specific application 
to Omega Gen. Either way, there do appear to be two 
important differences compared to our analysis in Sec- 
tion [121 First, NGB08 applied no corrections for incom- 
pleteness. We believe that this is generally inadvisable 
when using star counts, although it is not clear whether 
this may have specifically affected their analysis. Second, 
the NGB08 center was measured from the same GO-9442 
data set that we used, but from the text of their paper it 
appears that they used only the central pointing of the 
mosaic shown in Figure [1^. This pointing was centered 
on the traditional center of Omega Gen, which could con- 
ceivably have introduced some bias. More importantly, 
the use of only the central pointing implies that their ra- 
dial coverage went out to only 1!5, which is significantly 
less than the core radius. As a result, their analysis may 
have been more sensitive to the peak of the density distri- 
bution rather than to the symmetry point. This appears 
supported by the fact that their method did in fact iden- 



tify a density enhancement, as illustrated by Figure [20l 
This prompts the question whether there is any physical 
significance to this enhancement. 

The central part of Omega Gen has a large core that 
is almost homogeneous. Within this core, spatial gradi- 
ents in quantities of interest tend to very shallow. When 
only a finite number of stars are present or observed, one 
expects shot noise to dominate the small-spatial density 
distribution within the core. Some areas may be under- 
dense while others may be over-dense. This is in fact 
what we appear to be seeing. At the position that we 
have identified as the symmetry point of the cluster, the 
density appears somewhat underdense compared to its 
surroundings (see Figure [TO)) . By contrast, at the po- 
sition that NGB08 identified as the cluster center, the 
density appears somewhat overdense compared to its sur- 
roundings (see Figure [201) ■ However, these features have 
the hallmarks of being noise-induced valleys and peaks. 
As such, their presence has no bearing on where the ac- 
tual cluster center is. 

The apparent underdensity at our center, as well as 
the apparent overdensity at the NGB08 center, appear 
most prominent when the shot noise is highest. In Fig- 
ure \Wh this is because the number of stars in the mag- 
nitude range -6435 < —13 is low. In Figure [20b it is 
because of the use of unresolved light (which emphasizes 
the shot noise from the brightest stars, which are small in 
number). By contrast, in Figures [T^.c and Figure [^Ob 
the respective under- and over-densities are less promi- 
nent. Moreover, the errorbars are such that the profiles 
are consistent with being flat to within the uncertain- 
ties. Specifically, the rise in the number density within 
10" from the NGB08 center (Figure [20b) is not statisti- 
cally significant, and consistent with being a statistical 
fluctuation. The error bar on the central point is ~ 20%, 
and the point just outside the center is lower than the 
average within the inner 10". 

Even a simple visual inspection of Figure [TBI shows the 
core of Omega Gen is sparse and homogeneous enough 
to make identification of the center from local density 
enhancements either difficult or impossible. Instead, one 
should adopt a technique that focuses on the symmetry 
point of the larger-scale distribution. The data used for 
such an analysis must extend far enough out to cover and 
use the region where the density, brightness, and stellar 
motions start dropping significantly. From such analyses 
we have found that the NGB08 center is definitely not 
the symmetry point of the cluster. Also, any light or den- 
sity enhancement/cusp that may exist near the NGB08 
center appears consistent with a statistical fluctuation in 
an otherwise (nearly) homogeneous core. On top of this, 
we show in Section 16.21 below that the proper motions 
near the NGB08 center do not in any way indicate that 
this position is special compared to its surroundings. We 
therefore conclude that the position identified by NGB08 
is not the cluster center. 

6. ANALYZING THE PROPER MOTIONS 

The best way to constrain the presence of an IMBH 
in a cluster is to observe its effect on the motions of 
stars. NGB08 used an IFU on Gemini to measure the 
dispersion in the radial velocities of the unresolved light, 
and inferred a distinct rise in the velocity dispersion at 
their center, as compared with a field that was 14" away 
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Fig. 21. — (Top) The total proper motion ^jj = \J + /^^ 

for each star plotted as a function of distance from the center. 
(Bottom) The percentiles of the distribution for the 100-star bins 
indicated in the top panel. 

(see geometry in Figure [H]). We have demonstrated that 
the center they used is hkely 12" off from the true center. 
Nonetheless, if their velocity measurements are accurate, 
then they still imply an interesting kinematical feature 
within the core. To test this, we will examine the proper 
motions about their center and about the center we de- 
rived. 

6.1. Proper Motions in the Central 10 arcsec 

In Figure[2T]we plot the total motions no for the high- 
quality sample of stars (i.e., those flagged "good") in the 
central ten arcseconds as a function of distance from the 
cluster center. The total two-dimensional motions and 
errors are defined here as 



(1) 



The radial axis is scaled as so that we will get a roughly 
even distribution of stars across the graph. The vertical 
lines divide the sample into groups of 100 stars. No star 
in the inner 10" is moving faster than 3 mas/yr. 

This plot shows that the distribution of motions at the 
very center (leftmost bin) looks very similar to that in 
the last bin at ^--^10". There are no more high-velocity 
stars at the very center than elsewhere in the distribu- 
tion, contrary to what we would expect if an IMBH were 
present (Drukier & Bailyn 2003). The bottom panel ex- 
amines the percentiles of the bin-by-bin distributions. In 
addition to seeing no particularly fast-moving stars at the 
center, we also see no indication of a rise in the dispersion 
at the center. Consistent with our cursory examination 
above, the percentiles in the central bin at 2 arcsec- 
onds are indistinguishable from those in the outermost 
bin at 10 arcseconds. 

Figure indicates that the proper motions of Omega 
Cen do not show an obvious kinematical signature of an 



IMBH. By contrast, a sufficiently massive IMBH would 
have induced increasing velocities towards the center 
with RMS cr (X The question of what exact IMBH 

mass would be required to produce an observable signa- 
ture is discussed in detail in Paper II. One important 
issue when addressing this question is that many of the 
stars in the centermost parts of the field are not truly 
close to the center but are merely projected there from 
somewhere between ~ ±1 core radius along the line of 
sight. For example, for a projected aperture of 3" ra- 
dius around the center, only between 1.3% and 5.8% of 
observed stars reside within 3" from the center in three 
dimensions. These numbers were calculated from num- 
ber density distribution models derived in Paper II (the 
"core" and "cusp" models, respectively). Only 43 stars 
in the high-quality subset of our proper motion catalog 
reside within 3" from the projected center. Hence, for 
an IMBH mass that produces a sphere of influence of 
order 3", at most a handful of fast-moving stars would 
have been expected (see the quantitative analysis in Sec- 
tion 6.7 of Paper II). The bigger the IMBH mass, the 
more fast-moving stars would have been expected. 

6.2. Comparison with the NGB08 Kinematics 

The analysis in the previous section examined the mo- 
tions in the central region out to 10 arcseconds from our 
center, but it did not include the central field studied in 
NGB08. In Figure!^ we study the distribution of proper 
motions in the near vicinity of our center and the NGB08 
center, in the context of the stars in the wider central re- 
gion. The top row of plots shows the proper motions for 
the stars in the inner 15"; the middle and bottom rows 
show the same quantities for the stars within 3" of our 
center and the NGC08 centers, respectively. There are 
1200 stars in the wider central region, 43 stars near our 
center, and 54 stars near the NGB08 center. We include 
here only stars from the high-quality proper-motion sub- 
set. 

For each row, the leftmost panels show the proper- 
motion distribution as a vector-point diagram, and the 
right panels show the cumulative distributions. For the 
bottom two rows, we compare the cumulative distribu- 
tion for the given small region against the cumulative 
distribution for the central region as a whole. The dot- 
ted lines in the right plots show the result of 1000 Monte 
Carlo tests based on the distribution function from the 
upper panel and the observed number of stars in the 
lower two panels. The dotted lines indicate the 80% con- 
fidence region under the assumption that the velocity 
distributions are the same. 

We find that the motions of stars about our center and 
the NGB08 center are well described by the same distri- 
bution function that describes the stars within the inner 
15" of the cluster. Therefore, our proper motions do not 
confirm the velocity gradient reported by NGB08 in line- 
of-sight velocities, independent of the adopted center. 

To compare our kinematics more directly to those of 
NGB08, we also studied the proper motions of stars in 
exactly the same two 5" x 5" fields that they studied 
(see Figure [TB)) . In our high quality proper-motion sam- 
ple, there are 51 and 35 stars in their "central" and "off- 
center" fields respectively. Each star has two orthogonal 
velocity components, so the number of datapoints avail- 
able to estimate the one-dimensional velocity dispersion 
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Fig. 22. — A comparison of the proper-motion distribution within 
the inner 15" (top row of panels) and the distribution within 3" 
of the two centers (bottom rows). The left panels show a vector- 
point diagram, the middle panels show the cumulative distribution 
of the two dimensional proper motions fi£), and the right panels 
show the outer 25% of the cumulative distribution. The dotted 
lines in the bottom right plot indicate the 80% confidence region 
under the assumption that the underlying velocity distribution is 
the same as for the R < 15" region in the top panel. The observa- 
tions fall within the bands, indicating that there are no statistically 
significant differences between the distributions. 

cri_£) is twice the number of stars. Using the methods 
described in Paper II, we find that (Ji-d — 0.833 ±0.059 
mas/yr and 0.835 ±0.071 mas/yr, for the central and off- 
center fields, respectively. For a canonical Omega Cen 
distance D = 4.8 kpc (van de Ven et al. 2006), these re- 
sults translate to 18.9 ± 1.3 km/s and 19.0 ± 1.6 km/s, 
respectively. The similarity between these dispersions is 
consistent with the fact that we have determined both 
to be at roughly the same distance R — 12" from the 
center (see Figure [T6|) . These proper-motion dispersions 
are consistent with the values we measure throughout 
the rest of the central 15" region of Omega Cen (see 
Paper II). NGB08 determined line-of-sight velocity dis- 
persions (Tios = 23.0 ± 2.0 km/s and a\os = 18.6 ± 1.6 
km/s for these fields, respectively. Whereas our results 
agree with their measurement for the off-center field, our 
proper motions do not confirm the NGB08 result that 
the velocities in central field are higher. 

6.3. Radial Proper-Motion Profile 

Figure [23] shows the proper- motion distribution for the 
entire central data set, defined by the overlap region be- 
tween the GO-9442 and GO-10775 data sets (see Fig.H]), 
focusing on the percentiles of the distribution function 
(as marked on the plot). The plot on the left is more 
finely sampled and goes out to a radius of 50 arcseconds. 
The plot on the right distills the stars into 2000-star bins 
and goes out to the corners of the central field, at 120". 
The radial axis is once again scaled by , so that we will 
get equal number of bins with radius (where there is full 
azimuthal coverage) . 

While there is a distinct rise of about 10% in the PM 
dispersion from 80" into the center, the left plot shows 
that there is very little increase within the inner 30". 
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Fig. 23. — (Left) Similar to the bottom of Figure |2T] but cov- 
ering a larger radial extent. We take the stars 500 at a time, and 
determine the percentiles of the proper-motion distribution, which 
we plot as a function of the median radius in the group. (Right) 
Same, for 2000-star bins, out to the edge of the inner data set. In 
a two-dimensional distribution, 50% of the points should be within 
1.177 times the one-dimensional sigma (cri^jj), and 75% should be 
within 1.665 ai—o- 

This is as true for the wings of the distribution func- 
tion (the black symbols) as it is for the core (the cyan 
symbols) . 

In Paper II we will present a detailed analysis of 
the proper motions as a function of radius, including 
a determination of the velocity-dispersion and velocity- 
anisotropy profiles, comparison to literature data, calcu- 
lation of higher-order Gauss-Hermite moments, a com- 
parison between the proper motions in our central and 
major-axis fields, and a comparison between major-axis 
and minor-axis proper motions. Dynamical models will 
be fit to the data to constrain the possible presence and 
mass of any IMBH. 

7. EXAMINING THE PROPER MOTIONS BY 
POPULATION 

The primary motivation for the construction of a 
proper-motion catalog was to enable us to evaluate the 
likelihood of an intermediate-mass black hole at the clus- 
ter center (which we will further discuss in Paper II). 
However, this rich data set also can tell us much about 
the dynamical state of the cluster. In this section, we 
will do a cursory analysis of the data to examine the mo- 
tions of stars with different masses and stars in different 
populations. 

7.1. Equipartition 

It is well known that Omega Cen has not had enough 
time for complete dynamical relaxation. Harris (1996) 
reports its half-mass relaxation time as 10^° years. An- 
derson (2002) examined the luminosity function at the 
cluster center and at a removed radius and found that 
the cluster does not demonstrate as much mass segrega- 
tion as one would expect for a multi-mass King model in 
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PM IN EACH COORDINATE (mnn/yr) 

Fig. 24. — (Left) CMD showing stars selected for main-sequence 
study (those within the drawn-in hues). (Middle) The two- 
dimensional proper-motion errors cr^^ for stars as a function of 
instrumental F435W magnitude. Stars left of the drawn-in curve 
were used in the analysis. (Right) The x and y proper motions for 
the selected stars plotted together. The curves indicate the RMS 
as determined from a maximum-likelihood analysis that took into 
account the measured errors. The single-component RMS for each 
half-magnitude bin is reported on the right. 

energy equipartition. To complicate the matter, Omega 
Cen has multiple populations of stars, and these popu- 
lations exhibit spatial gradients and could have different 
mass functions. The present data set allows us to study 
the dynamics of the stars at a single location in the clus- 
ter (the center), thus avoiding the ambiguity introduced 
by comparison of populations at different radii. 

In Figure I24( we show the distribution of the proper 
motions and proper-motion errors as a function of 
F435W magnitude. For this analysis, we used the full 
PM catalog from the central field, not just the high- 
quality subset, since we wanted to examine the motions 
of the faint stars. The color-magnitude diagram on the 
left highlights the stars that we selected for this study 
of the main sequence (MS). The middle panel shows the 
distribution of PM errors for the stars that lie within the 
curves in the left panel. We selected for further analy- 
sis the stars that follow the general trend of PM error 
with magnitude in the middle panel, and show the mo- 
tions for these stars on the right. Both x and y motions 
are shown in the same plot. Within each half-magnitude 
bin, we found the error-corrected one-dimensional RMS, 
given the observed motions and errors for the stars in the 
bin. We report the RMS on the right-hand side of the 
plot. The proper-motion dispersion clearly increases as 
we go down the main sequence. The RMS motions dis- 
cussed here, as well as other kinematical quantities dis- 
cussed in the remainder of the paper, were determined 
using the maximum-likelihood methodology described in 
Appendix A of Paper II. 

In order to interpret this velocity variation in terms of 
mass, wc fit the CMD with an isochrone in the left panel 
of Figure 1251 finding a reasonable fit to the upper pop- 
ulation with a 12.5 Gyr isochrone from Pietrinferni et 
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Fig. 25.— (Left) The fit of our CMD with the [Fo/H] = -1.5 
and Y = 0.24 isochrone from Pietrinfirni et al. (2006). (Middle) 
The mass-luminosity relation for the isochrone, with symbols at 
the locations of our luminosity-function bins. (Right) The ob- 
served error-corrected one-dimensional proper-motion dispersion as 
a function of mass. The lines are drawn for comparison to indi- 
cate the trend that would be expected for complete equipartition 
and for no equipartition at all. The dotted line shows 
the intermediate case of a oc M~^'^. Solid points show measure- 
ments for the central field; open points show measurements for the 
adjacent major-axis field, renormalized upward by a factor 1.22. 



al. (2006) with alpha-enhancement, [Fe/H] = —1.6 and 
Y=0.24. This isochrone allows us to associate a mass 
with each F435W magnitude. The masses for our mag- 
nitude bins are shown in Panel (b). Finally, in Panel (c) 
we show (solid points) the run of error-corrected RMS 
proper motion with mass. If the cluster is in energy 
equipartition, we would expect the observed points to 
follow the upper curve, where velocity is proportional 
to 1/\/M. We see that the velocity does rise with de- 
creasing mass, but not as rapidly as equipartition would 
predict. These results are consistent with the core being 
in the process of establishing equipartition, but only be- 
ing about half-way there. Of course, it is well known that 
Omega Cen cannot be represented by a single isochrone. 
Combined with uncertainties in the exact distance, age, 
and metallicity, this causes small uncertainties in the 
absolute calibration of the mass for a given luminos- 
ity. However, this does not affect the conclusions about 
equipartition, which depend only on relative masses. 

We have repeated the central-field equipartition anal- 
ysis also for the adjacent major axis field. The median 
radius for stars in the central-field catalog is 75", whereas 
it is 218" for the major-axis field catalog. As a result of 
this difference in distance from the cluster center, wc find 
the RMS proper motion to be higher by a factor of 1.22 in 
the central field than in the major-axis field. However, 
Figure [^5b shows that the dependence of RMS proper 
motion on stellar mass is the same in both fields (open 
points indicate the error-corrected RMS proper motion 
measurements for the major axis-field after renormaliza- 
tion by a factor 1.22). Therefore, we detect no signifi- 
cant difference in the amount of equipartition between 
approximately 0.5 and 1.5 core radii. In principle, one 
might have expected less equipartition at the larger ra- 
dius, given that the two-body relaxation time increases 
with radius in Omega Cen (see figure 21 of van de Ven 
et al. 2006). On the other hand, the radial range that 
we can probe is not large. Also, cluster rotation be- 
comes a significant factor in the outer field, and this is 
not probed by our analysis (as discussed in Section [3.6.4p . 
Detailed evolutionary models would therefore be required 
for quantitative interpretation of these results, which is 
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outside of the scope of the present study. 

In Paper II we will compare the observed dispersions 
of line-of-sight velocities and proper motions in Omega 
Cen. The former are measured in km/s, while the lat- 
ter are measured in mas/yr. With the aid of dynam- 
ical models, this provides a means of determining the 
cluster distance. However, mass segregation provides an 
added complication in such an analysis. The stars for 
which we have proper motions are slightly less massive 
than the stars for which we have radial velocities (al- 
most all are giant stars). So even in the idealized case 
of an isotropic system, one would not expect to measure 
transverse and line-of-sight velocity dispersions that are 
the same in km/s. Figure [Ml allows us to determine the 
size of this effect and correct for it, as discussed in Paper 
II. 

7.2. Bulk motion of the metal-rich population 

Pancino et al. (2000) observed the RGB of Omega Cen 
from the ground in B and / and found evidence for three 
sub-populations: a metal-poor population (RGB-MP), 
an intcrmediatc-mctallicity population (RGB-Mint) and 
a metal-rich population (RGB-a), which comprises 5% of 
the RGB stars. Ferraro et al. (2002) then cross-identified 
these "RGB-a" stars in the proper-motion catalog of van 
Leeuwen et al. (2000) and found that they appeared to 
be moving at about 1 mas/yr relative to the rest of the 
cluster. They conjectured that this could represent a 
background cluster that may be in the process of merging 
with Omega Cen. This interpretation was disputed by 
Platais et al. (2003), who contended that it was likely to 
be a consequence of an uncorrected color-magnitude term 
in the plate equation. Recently, Bellini et al. (2009a) 
have constructed proper motions for stars in the outer 
regions of the cluster from ground-based CCD data span- 
ning four years and determined that the RGB-a stars do 
in fact share the bulk motion of the cluster. We will 
use our HST-measured proper motions and population- 
identifications to provide an additional determination of 
whether the metal-rich population is moving with the 
cluster. We focus here on the lower-turnoff stars, the 
SGB analog of the RGB-a stars. 

The left panel of Figure [25] shows a close up of turnoff 
region for the stars for which we have good proper mo- 
tions. We highlight the lower-turnoff (LTO) popula- 
tion in blue, and a control sample with similar flux in 
-Rf435W (and thus similar astrometric errors) in green. 
On the right, we show the proper motions for the two 
samples. It is clear that both distributions are centered 
on zero, meaning that they both share the bulk motion 
of the cluster. The cyan dot in the upper right panel 
shows the motion determined by Ferraro et al. (2002). 
The center of the distribution for the LTO popu- 
lation is —0.046 ± 0.040 mas/yr, and the fiy center is 
—0.062 ± 0.040 mas/yr. Both are consistent with zero, 
and similarly consistent with the motion for the control 
population. We therefore conclude that all of the popu- 
lations in the CMD are moving with the cluster and are 
phase-mixed. We note that the dispersions of the two 
populations are also the same, to within measurement 
errors (3%). 

7.3. Motions for the different MS populations 
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Fig. 26. — (Left) Color-magnitude diagram identifying the lower- 
turnoff population (blue in the on-line version) and the control 
population (red), which should have similar proper- motion errors. 
(Lower right) The vector-point diagram for the control sample. 
(Upper right) The vector-point diagram for the LTO population. 
The circle represents the 90%-confidence region for our fit to the 
center of the distribution. The filled dot (cyan) shows the motion 
as measured by Ferraro et al. (2002). 

Omega Cen was the first of the traditional globular 
clusters found to have multiple populations. The spread 
in metallicity along the giant branch has been known for 
many decades (Dickens & WooUey 1967 and Freeman & 
Rodgers 1975), but only recently have we been able to 
trace the multiple populations down to the unevolved 
stars. Anderson (2002) found and Bedin et al. (2004) 
confirmed that the main sequence clearly bifurcates into 
a red and blue branch below V ~ 20. Unexpectedly, 
they found that the fractions of stars in the two popu- 
lations are opposite to what one would expect based on 
the RGB populations and standard enrichment scenarios, 
which would have the bluer main-sequence population 
(bMS) being more metal poor (and more populous) than 
the redder population (the rMS). Norris (2004) pointed 
out that this could be explained by assuming that the 
metal-richer population was super-enriched in Helium. 
Villanova et al. (2007) then measured metallicities for 
stars in the two populations and found that the bluer 
stars were indeed more metal rich than the redder stars, 
in line with the Helium explanation. 

Since these initial discoveries of the split MS popula- 
tions, the spatial distribution of the stars has been stud- 
ied by SoUima et al. (2007) and by Bellini et al. (2009b). 
They find that the intermediate-metallicity population 
is more concentrated than the metal-poor population, 
both when the populations are identified on the RGB 
and on the main sequence. Kinematically, it was ini- 
tially believed that the M-Int population did not share in 
the cluster rotation (Norris & Freeman 1997), but recent 
spectroscopy by Pancino et al. (2007) mentioned above 
finds that all three populations appear to share the same 
rotation, to within the 2 km/s measurement errors. 

The radial velocities in the studies cited above come 
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Fig. 27. — Dispersion for the main-sequence stars the central 
field. (Left) CMD before selection, (Middle) CMD with selected 
bMS and rMS highlighted. (Upper right) PM distribution of the 
bMS stars, with the dispersions indicated at the top of the plot. 
(Lower right) Same for the rMS stars. Dashed circles are drawn 
for reference. The error in each dispersion is about 0.006 mas/yr. 

from all over the cluster, so they are able to probe the 
motions of stars with a global perspective. Here, we have 
proper motions in only two fields, one at the cluster cen- 
ter and one at about a 1.5 core radius out, along the 
major axis. Furthermore, we have motions only for stars 
on or below the SGB, so it is not possible to directly com- 
pare our motions against the radial velocities. Neverthe- 
less, we can still compare the motions of the populations 
we have access to. In Figures [27l and \28\ we compare the 
proper motions of the bMS and rMS stars in the central 
field and the major-axis field, respectively. 

Figure [27] shows motions for the bMS and rMS stars 
in the central field. The left panel shows the CMD in 
the location where the two MS populations are clearly 
distinguishable in the central field, between an instru- 
mental F435W magnitude of -11.4 and -10.6 (S/N - 
250). The middle panel shows our selection of the bMS 
and rMS stars. Finally, on the right, we show the vector- 
point diagrams with the motions for the two popula- 
tions. We analyzed the two distributions to determine 
the (error-corrected) dispersions as marked at the top of 
each plot. The dispersions for the two populations in 
both the major- and minor-axis directions are the same 
to within the measurement errors (0.006 mas/yr). For 
both MS populations, the major-axis dispersion is about 
3% greater than the minor-axis dispersion. All of the 
distributions have no mean motion, to within the mea- 
surement errors. 

Figure [28] shows a similar plot for the adjacent field, 
which is centered at a radius of about 4' (^--^1.5 Tc) along 
the SE major axis (see Fig[T]). Again, we identified the 
bMS and rMS stars at a brightness along the MS where 
we could clearly distinguish them and, at the same time 
where we had good proper motions. The third column of 
panels shows the vector-point diagrams for the two pop- 
ulations. Again, the error-corrected dispersions noted at 
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Fig. 28. — Similar to Figure (27] but for the adjacent field. The 
error in each listed dispersion is ~0.007 mas/yr. 



the top show no statistically significant differences be- 
tween the populations. 

The dispersion on the major axis is larger than that 
on the minor axis. This is what might naively have been 
expected from the fact that the dispersion provides the 
pressure that supports the shape of the system (rotation 
may generally contribute pressure as well, but we show 
in Section 5.2.3 of Paper II this is negligible near the 
center of the cluster). However, detailed axisymmetric 
anisotropic modeling would be required to fully inter- 
pret the differences between major axis and minor axis 
motion. 

It bears repeating that any mean motion was re- 
moved from the proper motions during the data reduc- 
tion stages, as discussed in Section 13.6.41 Hence, the 
motions as we have measured them here are unable to 
measure rotation directly. For that, one would need ref- 
erence objects to measure against (e.g., Anderson & King 
2003 measured the rotation of 47 Tuc in the plane of the 
sky using the background SMC stars). However, if there 
are multiple populations within the cluster and one pop- 
ulation happens to be rotating relative to another, then 
we would expect this relative rotation to manifest itself as 
a bulk motion between the two populations. This would 
be most prominent in the major-axis field, because ro- 
tation in stellar systems typically decreases towards the 
center. Of course there would be no rotation in the plane 
of the sky if the cluster were edge-on, but van de Ven et 
al. (2006) found the inclination of Omega Cen to be 48°. 
Therefore we would expect some difference in observed 
bulk motion if different populations had different rota- 
tions. 

Our analysis shows that the centers of the bMS and 
rMS distributions are the same in the major-axis di- 
rection, but differ by 3-a (0.037 mas/yr, or 0.8 km/s, 
assuming a distance of 4.8 kpc) in the minor-axis di- 
rection. This is consistent with the possibility that the 
populations are rotating relative to each other, but the 
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statistical significance of the result is only marginal (es- 
pecially when taking into account the possibility of small 
residual systematic effects at levels below 1 km/s). We 
similarly examined the lower-turnoff population in the 
major-axis field and found that it shows no systematic 
motion relative to the other populations. These findings 
are in agreement with the recent radial- velocity study by 
Pancino et al. (2007), who found that all the stars have 
the same rotational properties to within 2 km/s. 

The motions we have measured here provide only an 
incomplete picture into the dynamical state of the clus- 
ter. Ideally, we would like to have both accurate proper 
motions and radial velocities, as well as population- 
identification information, for a large number of stars 
distributed throughout the cluster. This would allow 
us to distinguish what kinds of orbits are populated by 
the different stellar populations, giving us much more 
information than simple studies of spatial distributions, 
dispersions, or rotational components. To this end, we 
could augment our proper-motion catalog with motions 
for the brighter stars using the shorter exposures; there 
exist measured radial velocities for many of these stars 
and we would then have 5 of the 6 phase-space coordi- 
nates (lacking only the line-of-sight location). However, 
the size of the central field would still be a significant 
limitation to our ability to infer the global properties of 
the cluster. 

8. CONCLUSIONS 

We have performed a careful reduction of the large 
GO-9442 data set which mosaic-imaged the inner 10' x 
10' of Omega Cen with well-dithered HST/ ACS obser- 
vations. We constructed a 1.2-million-star photometric 
catalog of positions and -Bf435w and i?F625W magnitudes 
from these data, along with a 14000 x 14000-pixel stacked 
image of the field in the same reference frame as the cata- 
log. We also reduced the data from two other ACS/WFC 
programs that overlap with the GO-9442 mosaic, but 
which covered only a single pointing (with medium-sized 
dithers). GO-10775 observed a central field and GO- 
10252 an adjacent field to the southeast, roughly along 
the major axis. Accurate astrometric analysis of the data 
allowed us to determine proper motions for stars in these 
fields. The resulting proper-motion catalogs, with cuts 
applied to retain only the "high quality" measurements, 
contain 53,382 stars in the central field and 19,593 stars 
in the major-axis field. The data products from our study 
are made publicly available as part of this paper (see Sec- 
tion ET]). 

We analyzed the positions of the observed stars to de- 
termine the cluster center. For this we used two sep- 
arate methods, one based on isodensity-contour fitting, 
and one based on the so-called "pie-slice" method. In 
the latter method we took particular care to model the 
effects of incompleteness, and to correct for them. The 
cluster centers thus determined each have an error bar of 
~ 1", and they agree to within the errors. Upon use of 
stars in common with the 2MASS catalog to calibrate to 
absolute coordinates, we find the center to be at (a, S) 
= (13:26:47.24, -47:28:46.45). We also used our proper- 
motion catalog to determine the kinematical center of 
the cluster, defined as the symmetry point of the proper- 
motion dispersion field on the projected plane of the sky. 
This represents the first time that the kinematic cen- 



ter of a globular cluster has been accurately determined. 
Again, we used methods based on contours and pie-slices, 
with consistent results. The kinematical center agrees 
with the star count center to within its 2" uncertain- 
ties. And finally, we also determined the center of unre- 
solved light in 2MASS data, again yielding a consistent 
result at the ~ 2" level. 

We computed the surface number-density profile of the 
cluster around its (newly determined) center. Artificial- 
star tests were used to correct for the effects of photomet- 
ric incompleteness. Density profiles were determined for 
various ranges of stellar magnitude, but were generally 
found to be similar, independent of magnitude. A single- 
mass King model provides a reasonable fit to the in- 
ferred profiles. There is no evidence for a strong number- 
density cusp towards the center; in fact the density in 
the centermost bins appears smaller (at marginal signif- 
icance) than that at somewhat larger radii {R « 10"). 

The proper-motion dispersion increases gently inwards 
from the core radius (~ 2.5 arcmin) to about 30", but 
flattens out at smaller radii. Detailed analysis in the cen- 
tral 15" shows little variation in kinematics with position. 
The dispersion does not increase appreciably towards the 
center, and the wings of the proper-motion distribution 
do not become more extended towards the center. There 
are no high-velocity outlier stars near the center that 
might be indicative of motion around an IMBH. 

We examined the variation of velocity dispersion with 
mass along the main sequence and found that although 
the dispersion does increase for the lighter stars, the clus- 
ter is not yet in equipartition. This is in agreement with 
the findings in Anderson (2002), who found that the clus- 
ter does not exhibit the mass segregation that would be 
expected for a multi-mass King model in equipartition. 
These results are not surprising, given the long half-mass 
relaxation time of ^ ^.O^'^^^'' ^'^ years (McLaughlin & van 
der Marel 2005). 

Omega Cen has long been known to have multiple stel- 
lar populations. The proper- motion catalog we have con- 
structed here has enabled us to look for variations in 
kinematics between populations. The blue (bMS) and 
red (rMS) main sequence stars have very similar dynam- 
ical properties. This holds for both the central field and 
the adjacent field along the major axis, and is true both 
in terms of dispersion and mean rotation (the latter is 
not a trivial result, because Omega Cen is not believed 
to be edge-on; van de Ven et al. 2006). The similarity in 
rotational properties contrasts with the findings of Norris 
& Freeman (1997) based on radial velocities. However, it 
agrees with a more recent study by Pancino et al. (2007). 
The mean motion of the metal-rich turn-off population 
is consistent with that of the rest of the cluster. This is 
what would be expected for a (quasi-)equilibrium config- 
uration, and disagrees with an earlier finding of Ferraro 
et al. (2002). 

NGB08 argued for the presence of an IMBH in the 
center of Omega Cen based on a combination of two 
arguments. First, they measured the integrated line-of- 
sight velocity dispersion of unresolved light in two 5" x 5" 
fields, one believed to be on the cluster center and one 
at i? = 14" from the center. The dispersion in the cen- 
tral field (23.0±2.0 km/s) exceeded that in the off-center 
field (18.6±1.6 km/s), consistent with the presence of an 
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IMBH of mass 4.0t°;o^ ^ 10^ Mq. Second, they measured 
the surface brightness profile of unresolved light. They 
found it to have a shallow central cusp of logarithmic 
slope 7 = 0.08 ± 0.03, consistent with theoretical predic- 
tions for the cusp induced by an IMBH (Baumgardt et 
al. 2005). The results from our new study have allowed 
us to test these arguments. 

We determined the one-dimensional proper- motion dis- 
persion of the stars in our catalog in each of the two 
fields studied by NGB08. For a canonical Omega Cen 
distance D = 4.8 kpc (van de Ven et al. 2006), the re- 
sults translate to 18.9 ± 1.3 km/s and 19.0 ± 1.6 km/s, 
for the central and off-center fields, respectively. So we 
find no kinematical difference between the fields, and we 
also do not detect kinematical gradients elsewhere in the 
central ~ 15". Moreover, we find that NGB08 did not 
actually observe the cluster center in their integral-field 
spectroscopy. The cluster center identified by NGB08 is 
12" from our newly determined center. We demonstrate 
that this (and similar offsets in other previous determina- 
tions) is likely due to biases induced by over-weighting of 
the small number of bright giants, or the limited region 
over which stars were measured (the central ACS chip 
in the case of NGB08, which covers only half the core). 
Here we have included the more plentiful main sequence 
stars and we have measured them out to beyond one core 
radius. This yields a precise and unbiased handle on the 
center of the cluster density distribution, as confirmed by 
our independent center determination from HST proper 
motion and 2MASS unresolved light data. 

The existence of a density cusp in Omega Cen, as 
reported by NGB08, has also not been confirmed by 
our analysis. Since they calculated their density profile 
around an incorrect estimate of the center, it is unclear 
what physical meaning their result may have. Either 
way, we showed here that use of unresolved integrated 
light, as was done by NGB08, is not the best way to 
constrain the cluster density profile. This method does 



not primarily measure the fiux from a large number of 
unresolved stars, but instead is sensitive to the large- 
radii scattered PSF wings of bright giants. This method 
therefore suffers more from shot noise (biases related 
to small-number statistics) than a measurement of the 
number-density profile. This is particularly important 
in a (nearly) homogeneous core such as that in Omega 
Cen, since shot noise will always cause some areas of high 
apparent density to exist by chance. The error bars on 
our number density profile around the NGB08 center are 
such that the apparent density enhancement there is con- 
sistent with being a statistical fluctuation. The number 
density profile determined by us around the newly de- 
termined center also shows little evidence of a significant 
cusp towards the center, although a shallow cusp may 
not be ruled out. 

In summary, our results do not validate the arguments 
put forward by NGB08 to suspect the presence of an 
IMBH in Omega Cen. However, this does not mean that 
such an IMBH may not be present after all. Our new 
proper-motion catalog far exceeds the quality and quan- 
tity of the kinematical data in the central arcmin previ- 
ously available. This provides the opportunity to study 
the central dynamics of Omega Cen at a level of detail 
that is unmatched by almost all other clusters, with the 
possible exception of 47 Tuc (McLaughlin et al. 2006). 
In Paper II we therefore present a new detailed study of 
the dynamics and density profile of Omega Cen, with the 
primary goal of exploiting the new data to constrain the 
mass of any possible IMBH. 
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TABLE 6 

Various centers. (The NGB08 center has been corrected for HST guidestar errors.) 



Ref 


RA 


Dec 


X 


y 


Ax 


Ay 


Ar(") 


This work 


13:26:47.24 


-47:28:46.45 


6725 


6810 








NGB08 


13:26:46.08 


-47:28:42.9 


6962 


6881 


+237 


+71 


12.3 


Harris 


13:26:45.9 


-47:28:36.9 


7011 


6960 


+286 


+150 


16.1 


van Leeuwen 


13:26:45.756 


-47:28:42.78 


7027 


6883 


+302 


+73 


15.5 
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